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ABSTRACT 

Context, The [C n] 158yum line is an important tool for understanding the life cycle of interstellar matter. Ionized carbon is present in 
a variety of phases of the interstellar medium (ISM), including the diffuse ionized medium, warm and cold atomic clouds, clouds in 
transition from atomic to molecular, and dense and warm photon dominated regions (PDRs). 

Aims. Velocity-resolved observations of [C n] are the most powerful technique available to disentangle the emission produced by 
these components. These observations can also be used to trace CO-dark Hi gas and determine the total mass of the ISM. 
Methods. The Galactic Observations of Terahertz C+ (GOTC+) project surveys the [Cn] 158/um line over the entire Galactic disk 
with velocity-resolved observations using the HerschelfHIFl instrument. We present the first longitude-velocity maps of the [C n] 
emission for Galactic latitudes b = 0°, ±0.5°, and ±1.0°. We combine these maps with those of Hi, 12 CO, and 13 CO to separate the 
different phases of the ISM and study their properties and distribution in the Galactic plane. 

Results. [Cn] emission is mostly associated with spiral arms, mainly emerging from Galactocentric distances between 4 and lOkpc. 
It traces the envelopes of evolved clouds as well as clouds that are in the transition between atomic and molecular. We estimate that 
most of the observed [C n] emission is produced by dense photon dominated regions (~47%), with smaller contributions from CO- 
dark H 2 gas (~28%), cold atomic gas (~21%), and ionized gas (~4%). Atomic gas inside the Solar radius is mostly in the form of 
cold neutral medium (CNM), while the warm neutral medium (WNM) gas dominates the outer galaxy. The average fraction of CNM 
relative to total atomic gas is ~43%. We find that the warm and diffuse CO-dark H2 is distributed over a larger range of Galactocentric 
distances (4-11 kpc) than the cold and dense H 2 gas traced by 12 CO and 13 CO (4-8 kpc). The fraction of CO-dark H 2 to total H 2 
increases with Galactocentric distance, ranging from ~20% at 4 kpc to ~80% at 10 kpc. On average, CO-dark H 2 accounts for ~30% 
of the molecular mass of the Milky Way. When the CO-dark H 2 component is included, the radial distribution of the CO-to-H 2 
conversion factor is steeper than that when only molecular gas traced by CO is considered. Most of the observed [C 11] emission 
emerging from dense photon dominated regions is associated with modest far-ultraviolet fields in the range xo — 1 - 30. 

Key words. ISM: atoms — ISM: molecules — ISM: structure — Galaxy: structure 

1. Introduction the warm ionized medium, the warm and cold diffuse atomic 

medium, and warm and dense molecular gas. The [Cn] 158yum 
The transition from atomic to molecular clouds, the formation line is the main coolant of the diffuse ISM and therefore plays a 

of stars within high density regions, the radiative feedback from key m]s in the thermal instability that conve rts from warm and 



O ■ newly formed stars, and the disruption of molecular clouds and diffuse atomic douds to cold and dense douds As the [c n] in _ 

•n ; termination of star formation are of great astrophysical inter- tensity is sensitive to column densityi volume densityi and kinetic 

7H , est as they are processes governing the evolution of galaxies temperature, it can be used to determine the physical conditions 

£>. ■ in our Universe. In the Milky Way the evolution of the inter- fth e line-emitting gas 
• ■ i-h stellar medium (ISM) has traditionally been studied with radio 

/\ , continuum and Ha observations for the ionized component (e.g. The first sensitive all sky survey of line emission in the wave- 
H ■ lHaffner et al.ll2009l) . with obse rvations of the Hi 21 c m line for length range between lOOyurn and 1 cm was car ried out with 
- - .' the atomic gas component (e.g. lKalberla & Kerpl2009t) . and with the Fa r-Infrared Spectrometer (FIRAS) on COBE dBennett et al.l 
rotational transitions of the CO molecule f or the dense molecu- Il994h . The angular resolution of COBE was 7° and the spec- 
lar gas component (e.g. iDame et al.ll200ll) . But observations of tral resolution of FIRAS was about 1000 kms 1 for the [Cn] 
the transition between these different stages of evolution with 158yum and [Nn] 205//m lines. The COBE map showed that 
velocity-resolved observations required to separate different gas [Cn] is widespread over the Galactic plane and that it is the 
components along the line-of-sight (LOS) have been missing. brightest far-infrared line. COBE also found a correlation be- 
The [Cn] 158/im fine-structure transition is an excellent tween the [Cn], [Nn], and far-infrared (FIR) continuum emis- 
tracer of the different stages of evolution of the ISM. As ionized sion, suggesting that [Cn] might be a good tracer of star- 
carbon (C + ) can be found throughout the ISM and can be excited form ation activity in galax ies. The balloon-borne BICE mis- 
by collisions with electrons, Hi, and H 2 , the [Cn] line traces sion dNakagawaet al.ll 1998b observed [Cn] in the inner Galactic 

plane with 15' angular resolution and 175 km s velocity resolu- 

* Herschel is an ESA space observatory with science instruments tion. This mission gave better constraints on the latitudinal distri- 

provided by European-led Principal Investigator consortia and with im- bution of the Galactic [C n] emission. The COBE and BICE mis- 

portant participation from NASA. sions did not have the required angular and velocity resolution to 
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resolve spatially and in velocity individual clouds in the Galactic 
plane. The Kuiper Airborne Observatory allowed the study of 
a handful o f Hn regions with velocity -resolved fC nj obser- 
vations (e.g IBoreiko et al.lll988t iBoreiko & Betzlll99ll) . These 
obse rvations, and also tho se recently obtain ed with Hersche l 
(e.g. lOssenkopf et al.ll20Tl and SOFIA (e.g. lGraf et aDl2012h . 
were mostly limited to hot and dense photon-dominated regions 
(PDRs) associated with massive star-forming regions. 

Far-infrared (FIR) dust continuum emission has been 
widely used as a mass tracer of the ISM, and several high- 
resolution, large-scale, sensitive maps o f the Galactic plane exist 
dCarev et al.ll2009t iMolinari et al.ll2010l) . Dust continuum emis- 
sion is the result of the integrated contribution of different ISM 
components distributed along the line-of-sight. Therefore, it is 
impossible to use this emission to isolate and study the different 
contributing ISM components, particularly for observations of 
the Galactic plane. 

The Herschel dPilbratt et al.ll2010l) Open Time Key Project 
Galac tic Observations of Terahertz C+ (GOTC+: lLanger et al.l 
l201ll) provides the first high-angular resolution (12"), velocity- 
resolved (0.1 km s" 1 ) observations of the [Cn] line in the 
Galactic disk, allowing us to resolve individual clouds and to 
separate the different ISM components along the line-of-sight. 
The GOTC+ survey consists of ~500 LOSs distributed in a 
volume-weighted sampling of the Galactic disk. We have com- 
plemented our GOTC+ [Cn] data with observations of the Hi 
21 cm line, tracing diffuse atomic gas, and of CO and its iso- 
topologues, tracing the dense molecular gas. 

In this paper we present the full GOTC+ Galactic plane sur- 
vey. We use the kinematic information provided by high ve- 
locity resolution observ ations of [Cn], obtained with the HIFI 
(Ide Graauw et al.lfeOlOl) instrument on Herschel, together with 
complementary CO and H 1 surveys to study the structure of the 
galaxy and to characterize the distribution of the different com- 
ponents of the interstellar medium in the Milky Way. 

Neutral atomic hydrogen represents the dominant gas com- 
ponent of the interstellar medium in galaxies. This H 1 gas is pre- 
dicted to_£xisHnJwopJiasesJnrough thermal pressu re equ ilib- 
rium (Pi kerNerll968l Field et al. 1969: IWolfire et all l995); the 
cold neutral medium (n H =* 50 cm" 3 ; T^ n ^80 K; CNM) and the 
warm neutral medium (« H - 0.5 cm 3 ; 7\ in ^8000 K; WNM). 
Note, however, th at gas in thermally unstable conditions has 
been observed (e.g lDickevetal.il 19771: iHeiles & Trolandll2003l) 
and p redicted in numerical simulations (see Va zquez-Semadenil 
120121 and references therein). The atomic gas component has 
been extensively mapped using the H 1 2 1 cm line in the Milky 
Way and external galaxies. But this line in emission traces only 
the total column of gas along the line-of-sight (assuming that it 
is optically thin), and therefore it is impossible to disentangle the 
relative contribution from gas in the CNM or WNM. H 1 absorp- 
tion studies of nearby clouds in the foreground of extragalac- 
tic continuum sources have provided constraints on the relative 
fraction of the CNM and WNM, sugge sting that 60% of the ga s 
is in the WNM and 40% is in the CNM dHeiles & Trolandl2003l) . 
But these observations are limited to a small number of sources 
distributed mostly at high Galactic latitudes, thus corresponding 
only to the local ISM. We therefore lack the knowledge about 
the distribution and conditions of the CNM and WNM over the 
entire Galactic disk. The [C 11] 158/im line is also a tracer of the 
diffuse neutral gas but has the advantage that it is also sensitive 
to the density and temperature of the gas. In particular, because 
of the density contrast between the CNM and WNM, for a given 
H 1 column density the [C 11] emission associated with the WNM 
is expected to be a factor of ~20 weaker than that associated 



with the CNM (IWolfire et al.ll2010h . We can therefore use the 
GOT C+ [C n] emission in combination with H 1 to study the dis- 
tribution and properties of the WNM and CNM over the entire 
Galactic disk. 

The mass of molecular clouds has traditionally been traced 
using observations of the CO molecule. However, the column 
density from the cloud exterior to the location where hydrogen 
becomes molecular is smaller than that required for carbon to 
become molecular in the form of CO. It is therefore expected 
that in the outer parts of molecular clouds a significant fraction 
of the H2 is associated with C + and C° instead of with CO. This 
gas is called "hidd en-H?" or "dark-H?" gas (see iGrenier et al.l 
120051: lLanger etal.ll201Ct IWolfire et al.ll2010t) . In this paper we 
will refer to this component as "CO-dark H2 ga s". Indirect ob- 
servations of ga mma-rays ( Qrenieret al. 2005) and dust con- 
tinuum emission dPlanck Collaboration et al.ll201 lb suggest that 
30-50% of the molecular mass in the Galaxy is not traced by 
CO. The preliminary results from GOT C+, based on a sample 
of 16 LOSs, reached a similar conclusion (Langer et al.lEOlOl 
IVelusamv et alJl2010l) . IVelusamv et al.l (1201 3l) used a Gaussian 
decomposition approach to study the GOTC+ LOS in the in- 
ner Galaxy, -90°< I < 57°, discussing the distribution of the 
[C n] velocity components as a function of Galactocentric dis- 
tance and the distribution of CO-dark H2 gas as function cloud 
type (diffuse clouds, transition clouds, and dense clouds traced 
by 13 CO). A more detailed analysis of the Gaussian decompo- 
sition of the GOT C+ [C n] emission in the inner galaxy will be 
presented in Langer et al. (2013) in preparation. Note that some 
of the H2 gas might not be traced by CO due to the lack of sen- 
sitivity of th e CO ob servations. In the Taurus molecular cloud, 
iPineda et all d2010al) found that 23% of the total mass is in the 
outskirts of the cloud where neither 12 CO nor 13 CO could be de- 
tected in single pixels, but both isotopologues could be readily 
measured by averaging pixels over a large area. 

A significant fraction of the FIR line and contin- 
uum emission in galaxies originates from dense photon- 
dominated regions (or phot odissociation regions; or PDRs; see 
iHollenbach & Tielensl 19991 and references therein). In PDRs the 
chemistry and thermal balance is dominated by the influence of 
far-ultraviolet photons from massive stars. Therefore, PDRs are 
the best sites for studying the radiative feedback from newly 
formed massive stars on their progenitor molecular gas, a pro- 
cess that plays a key role in the regulation of star-formation in 
galaxies. The [C 11] line is a very bright tracer of PDRs and, to- 
gether with fine-structure transitions of [0 1] and [C 1] and rota- 
tional transitions of CO, can be used to determine the physical 
conditions of the line-emitting regions. The FUV field inten- 
sity is a key parameter for PDRs and is closely related to the 
massive star-formation activity arou nd PDRs. T he FUV field 
is often measured in units of the Draine ( 19781) fielqj defined 
as the average FUV intensity in the local ISM. iKaufman et al] 
(1 19991) showed that the ratio of [C 11] to CO is a good tracer of 
the strength of the FUV field in the ga l axy. Us ing this line ratio 
in a sample of 16 LOSs. IPineda et aD d2010bl) found that most 
of the dense PDRs have FUV fields lower than 100 times the in- 
terstellar FUV field. With the complete GOTC+ survey, we can 
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The Draine field is isotropic 



(i.e. a given point is illuminated from 4ir steradians), while the surface 
of the clouds considered here are only illuminated from 2n steradians, 
therefore the rate of photoreactions at the cloud surface are half of what 
they would be with the Draine field. 
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use the [C n]/CO ratio to trace the strength of the FUV radiation 
field over the entire Galactic disk. 

This paper is organized as follows. In Section[2]we describe 
the GOT C+ observations and data reduction. We also describe 
the observations of ancillary CO and H i data. In Section [3~Tl we 
present position-velocity maps of the [C n] emission and relate 
them to the spiral structure of the Milky Way. In Section 13.21 
we present the distribution of the [Cn], Hi, and CO emissivi- 
ties as a function of Galactocentric distance. In Section [4] we 
combine the [Cn], Hi, and CO emission to identify the differ- 
ent ISM components and determine their contribution to the ob- 
served [C n] emission. We describe the distribution of the warm 
and cold neutral medium atomic gas components in Section lBTTl 
In Section [5T2l we study the radial distribution of the CO-dark 
H2 gas component and compare it to that of molecular gas traced 
by CO and I3 CO In Section|6]we present a determination of the 
FUV field distribution in the Galactic plane. We summarize our 
results in Section [7] 



2. Observations 

2.1. [Cn] observations 

We surveyed the Galactic plane in the [C 11] 2 P3/2 — > 2 Pi/2 fine 
structure line at 1900 .5469 GHz (rest frequency) with the HIFI 
dde Graauw et al.l2010l) instrume nt on-board the Herschel space 
observatory dPilbratt et al.ll2010l) . These observations are part of 
the Herschel Open Time Key Project Galactic Observations of 
Terahertz C+ (GOTC+), which is devoted to studying the [Cn] 
158 fim line over a wide range of environments in the Galactic 
plane, the Galactic center, and in the Outer Galaxy Warp. The 
GOTC+ Galactic plane survey consists of 452 LOSs distributed 
in an approximately volume-weighted sample of the Galactic 
plane. We present a selection of GOT C+ lines-of-sight in [C 11], 
Hi, and 12 CO emission in Figure Q] The entire data set is pre- 
sented as a set of position-velocity maps in Figure [2] and in 
Appendix [B] We sampled in Galactic longitude every 0.87° for 
|/| < 60°, every 1.3° for 30°< \l\ < 60°, and every 4.5° for 
60°< |/| < 90°. In the Outer Galaxy (|/| > 90°), we sampled 
in alternating steps varying from 4.5° to 13.5°. We sampled in 
Galactic latitude 0.5° and 1.0° above and below the plane at 
alternate positions in /. In 11 LOSs in the Outer Galaxy with 
|/| > 90°, we sampled b = +2.0° instead of +0.5°. 

The [Cn] 1.9 THz observations were carried out with 
the HIFI Band 7b receiver, which is employs Hot Electron 
Bolometer (HEB) mixers. The HEB bands in HIFI show promi- 
nent electrical standing waves that are produced between the 
HEB mixing element and the first low noise amplifier. The stand- 
ing wave shape is not a standard sinusoid and is difficult to re- 
mov e from the result i ng sp ectrum using standard fitting meth- 
ods dHiggins & Kooil r2009h . To remove these standing waves 
we used a procedure developed at the HIFI Instrument Control 
Center (ICC; Ian Avruch priv. comm.) which generates a li- 
brary of standing wave shapes from different observations and 
finds the best fitting one to correct the observed spectrum (see 
lHigginsl201ll for a detailed description of this method). We pro- 
cessed our data with this standing wave removal script in HIPE 
dOtt et al.l 120061) version 8. In HIPE we also removed residual 
standing waves by fitting a single sinusoidal function using the 
FitHIFIFringeO procedure. After all standing waves are re- 
moved, we exported our data to the CLASS9S0 data analysis soft- 
ware, which we used to combine different polarizations, fit poly- 
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Fig. 1. A selection of GOT C+ lines-of-sight observed in [C 11], 
Hi, and 12 CO emission. Note that the [Cn] emission systemat- 
ically extends closer to the tangent velocity compared with the 
12 CO line, where it follows the distribution of the Hi emission 
(see Section [3~TI ). 



nomial baselines (typically of order 3), and smooth the data in 
velocity. 

The angular resolution of the [C 11] observations is 12". We 
applied a main-beam efficiency of 0.72 to transform the data 
from an antenna tempe rature to a main-beam temperature scale 
dRoelfsemaetalJl2012l) . 

We used the data produced by the wide band spectrometer 
(WBS), which have a channel width of 12 MHz (0.16 km s" 1 at 
1.9 THz). We later smoothed the data to a channel width of 0.8 
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km s ' . For this channel width the average rms noise of our data 
iflO.lK. 

As we expected the [Cn] emission to be extended in the 
Galactic plane, we used the LoadChop with reference observing 
mode. In this observing mode an internal cold calibration source 
is used as a comparison load to correct short term changes in 
the instrument behavior. Since the optical path differs between 
source and internal reference, a residual standing wave struc- 
ture remains and is corrected with an observation of a sky refer- 
ence position. For a given Galactic longitude, LOSs with b — 0°, 
+0.5° and +1.0° share the same reference position at b — ±2.0°. 
In the 11 Outer Galaxy positions where we observed b - 0°, 
+ 1.0°, and +2.0°, we used the b - +2.0° as a reference for 
b = 0°, +1.0° and b = ±4.0° as a reference for b = ±2.0°. We re- 
covered the reference spectrum (OFF) using HIPE by subtracting 
from the ON-OFF spectrum a spectrum for which the reference 
subtraction step in the HIFI pipeline was suppressed. The result- 
ing spectra include a standing wave with period ~90MHz. We 
combined all observations of the reference spectrum in CLASS90 
and removed the standing wave using a Fast Fourier Transform 
(FFT) procedure. In the cases where emission in the reference 
position is present, we decomposed this emission into Gaussian 
components and corrected the ON-source spectra appropriately. 
A total of 5 1 % of the GOT C+ spectra were corrected for emis- 
sion in the reference position. 



3 For the typical [C n] FWHM line width of about 3 km s (Langer 
etal. 2013 in preparation), this sensitivity limit corresponds to l.lxlO -6 
erg s cnT 2 sr" 1 . The integrated intensity in units of Kkrns -1 can be 



2.2. CO observations 

To complement the GOTC+ data, we observed the / = 1 — > 
transitions of 12 CO, 13 CO, and C 18 with the ATNF Moprfl 
Telescope. We observed all GOTC+ positions towards the inner 
Galaxy between / = -175.5°and / = 56.8° with an angular res- 
olution of 33". Typical system temperatures were 600, 300, and 
250 K for 12 CO, 13 CO, and C 18 0, respectively. To convert from 
antenna to main-b eam temperature scale, we used a main-beam 
efficiency of 0.42 (lLadd et alj|2005l) . All lines were observed si- 
multaneously with the MOPS spectrometer in zoom mode. The 
spectra were smoothed in velocity to 0.8 km s for 12 CO and 
I3 CO and to 1.6 km s 1 for C lo O. The typical rms noise is 0.6 K 
for 12 CO and 0. 1 K for both 13 CO and C 18 0. We checked point- 
ing accuracy every 60 minutes using the closest and brightest 
SiO maser. For 102°< I < 141°, we used th e 12 CO J = 1 -> 
data from the FCRAO Outer Galaxy survey dHeveret al.lll998l) . 
The data have an angular resolution of 45", and a rms noise of 
about 0.6 K in a 0.8 km s -1 channel. 



2.3. Hi data 

We used supplementary Hi 21 cm dat a taken from the Souther n 
Galactic Plane Survey (SGPS; iMcClure-Gr iffiths et al. 2005), 
the Canadian Galactic Plane Survey (CGPS: iTavlor etail 2003 ). 



and the VLA Galactic Plane Survey (VGPS: IStil et al .1120061) 
These surveys combined cover the 253°< I < 67° portion of 
the Galactic plane. The SGPS (253°< I <358°) data have an an- 
gular resolution of 2', and a rms noise of 1.6 K per 0.8 km s 
channel. The CGPS covers 74°< I <147°, with an angular res- 
olution of 1', and a rms noise of 2K in a 0.82kms~ 1 channel. 
Finally, the coverage of the VGPS is 18°< / <67°, with an an- 
gular resolution of 1', and a rms noise of 2 K in a 0.82 km s 
channel. For the Galactic center we used the data presented by 



converted that in units of erg s 

10 5 /[erg s" 1 cnT 2 sr 1 ] JGoldsmith et al.l2012h . 



cm sr u sing /[Kkms ] = 1.43 x 



4 The Mopra radio telescope is part of the Australia Telescope which 
is funded by the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. 



Pineda, Langer, Velusamy, and Goldsmith: A Herschel [C n] Galactic Plane Survey 



iMcClure-Griffiths et al.l (1201 2l) which have parameters similar to 
those of the SGPS . We also used the GALFA-H i Arecibo survey 
dPeeketal.ll2011h to cover the Galactic longitude range between 
1 80° < / <2 1 2° . The GALFA-H i data have an angular resolution 
of 4' and a typical rms noise of 80 mK in a 1 km s channel. 



3. The Distribution of [Ch] in the Milky Way 

3. 1 . Position-Velocity Maps 

In Figure [2] we present the GOTC+ survey in a set of position- 
velocity maps for b — 0°. The corresponding position- velocity 
maps for ±0.5°, and +1.0° are presented in AppendixlBl The ob- 
served [C n] emission peaks at b — 0° and decreases for LOSs 
above and below the Galactic plane. Most of the [Cn] is con- 
centrated in the inner galaxy (-60°< / <60°). The position- 
velocity map is overlaid with projections of the Scutum-Crux, 
Sagittarius-Carina, Perseus, and Norma-Cygnus£| Milky Way 
spiral arms. We used the fits to the parameters determining 
the log arithmic spiral arms presented bv lSteiman-Cameron et al.l 
(1201 Oh and assumed a flat rotation curve the distance of the Sun 
to the Galactic center, R G - 8.5 kpc, and an orbital velocity of the 
sun with respect to the Galactic center, V = 220 km s -1 , which 
are the values recommended by the International Astronomical 
Union (IAU). We show the projected spiral arms in the range 
of Galactocentric distance between 3 kpc and 25 kpc. We can 
see that most of the [C n7 emission is closely associated with 
the spiral arms and it is brightest at the spiral arms' tangent 
points. In Figure [3] we plot the spiral arms with a color scale 
that represents the Galactocentric distance. Comparing Figure[3] 
with the panels in Figure [2] we see that the [C n] emission is 
mostly concentrated within the inner 10 kpc of the Galaxy (see 
also Section l3J2l . The [Cn] emission is slightly shifted with re- 
spect to that of 12 CO, with the [C n] emission present at veloci- 
ties that are closer to the tangent velocity. This shift can also be 
seen near the tangent velocities in the sample spectra shown in 
Figure Q] This suggests that the [C n] emission traces an inter- 
mediate phase between the extended H i emission and the cold 
and dense molecular gas traced by 12 CO. 



3.2. Radial Distribution Plots 

In the following we study the distribution of the [C n] emission 
as a function of Galactocentric radius. We divided the galaxy 
in a set of rings with a distance to the Galactic center R ga \ and 
width AR - 0.5 kpc. For each of these rings, we calculated 
the azimuthally-averaged integrated intensity per unit length or 
emissivity of the [C n] line. We describe the procedure used to 
determine these emissivities in AppendixlAl 

In Figure |4] we present the azimuthally averaged emissivity 
of [Cn] as a function of Galactocentric radius. The [Cn] emis- 
sion shows a main peak at 4.2 kpc and a secondary peak at ~ 6 
kpc. As mentioned above, the [C n] emission is concentrated in 
the inner 10 kpc of the Galaxy. When we separate the contri- 
bution from LOSs with / = - 180° and / = 180 - 360°, we 
can see that the peak at 4.2 kpc originates from LOSs in the 
1 = — 180° range while the peak at 6 kpc originates from LOSs 
in the I — 180 - 360° range. This difference in the emissivity dis- 
tribution between both sides of the Galaxy illustrates the mag- 
nitude of the intrinsic scatter within the rings, which is a result 
of the fact that the Galaxy is not axisymmetric. Because the aim 
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5 The Norma-C ygnus arm is also known as the Norma-3kpc arm 
fe.g. lValleeil2008h . 



Fig. 4. Radial distribution of the azimuthally averaged [Cn] 
emissivity at b = 0°. The shaded area represents the emissiv- 
ity calculated considering all sampled Galactic longitudes. The 
connected boxes represent the radial emissivity distribution for 
0°< / < 180°, while the connected circles represent that for 
180°< / < 360°. The typical 1 <x uncertainty in the [C n] emissiv- 
ity is 0.02 K km s" 1 kpc" 1 . 



of our analysis is to describe the average properties of the ISM 
within a Galactocentric ring, we only consider the uncertainties 
that describe how accurately this azimuthal average emissivity 
is determined. These uncertainties are generally very small, as 
they are the result of averaging over a large number of data val- 
ues. In many figures the error-bars are too small to be visible, 
especially when the emissivity vary over a large range, and we 
will quote them in the figure captions. Note that in most cases 
the uncertainties in the assumptions of our modelling are much 
larger than those indicated by the observational error-bars. 

In Figure[5]we show the contributions to the [C n] emissivity 
from the different spiral arm tangents. We consider LOSs in 12°- 
wide bins in Galactic longi tude center ed at the spiral arm tan- 
gent locations defined by Vallee (2008). Note that different spi- 
ral arms can be seen in this set of LOSs, but the peak of the [C n] 
emission corresponds to the tangent of the spiral arm. The peak 
at 4.2 kpc is associated with the Scutum-Crux tangent while the 
6 kpc peak is related to the combined emission from the Start 
of Perseus (peak at 4.7 kpc), Norma-Cygnus (peak at 5.7 kpc), 
and Crux-Scutum (peak at 6.7 kpc) tangencies. The Sagittarius- 
Carina (peak at 6.2 kpc) and Carina-Sagittarius (peak at 8.7 kpc) 
make a smaller contribution. Figure [5] shows that the brightest 
[C n] emission in the GOT C+ survey arises from the Scutum- 
Crux tangent at 4.2 kpc. 

In Figure|6]we show the radial distribution of the [C n] emis- 
sivity from LOSs with b = 0°, b = +0.5°, and b = +1.0°, and 
for all observed LOSs together. The [Cn] peak at 4.2 kpc from 
the Scutum-Crux arm is only seen for b — 0°, while the intensi- 
ties from b = +0.5°, b = +1.0° mostly contribute to the inten- 
sity at 5.2 kpc from the Galactic center. The shift in the peak of 
the radial distribution for \b\ >0°is related to the vertical struc- 



14 



Pineda, Langer, Velusamy, and Goldsmith: A Herschel [C n] Galactic Plane Survey 

7 



T 1 1 

Start of Perseus 
Norma-Cygnus 
Crux-Scutum 
Carina-Sagittarius 
Sagittarius-Carina 
Scrutum-Crux 




5 6 7 8 
Galactocentric Distance [kpc] 



Fig. 5. Radial distribution of the azimuthally averaged [C n] 
emissivity fo r GOTC+ LO Ss toward the tangents of spiral arms 
as defined bv lValled d2008l) . 
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Fig. 6. Radial distribution of the azimuthally averaged [C n] 
emissivity for LOSs at b = 0°, b = +0.5°, b = ±1.0°, and for all 
observed LOSs together. 



ture of the Galactic disk. A cloud with R gi \ = 5 kpc observed 
in a LOS with |fo|=l° would have a vertical distance from the 
Galactic plane of ~65 pc, assuming that it is in the near side of 
the Galaxy, while a cloud with R sa i = 4 kpc, in the same situa- 
tion, would have a vertical distance of ~100pc. A Gaussian fit 
to t he latitudinal [C nl dist ribution observed by BICE (see Fig. 
9 in iNakapawa et aljfl998h results in an angular FWHM of the 
Galactic disk of 1 .9°. Assuming that most of their observed [C n] 
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Fig. 7. Radial distribution of the azimuthally averaged [C n], H i, 
and 12 CO emissivities at b = 0°. The discontinuity in the Hi 
emissivity distribution at ~9 kpc is discussed in Section I5.ll 
Typical \cr uncertainties in the [Cn], Hi, and 12 CO emissivities 
(before scaling) are 0.02, 0.6, and 0.1 Kkms^'kpc -1 , 
tively. 



respec- 



emission is in the near side of the Galaxy at a heliocentric dis- 
tance of about 4 kpc, this angular FWHM corresponds to 130pc. 
Our inability to detect a peak of emission at 4 kpc for \b\ > 0° 
seems to be consistent with this value of the [C n] disk thick- 
ness. Results on the vertical structure of the Galactic disk using 
the GOTC+ data will be presented in a future paper. 

In Figure [7] we compare the radial distributions of the [C n], 
12 CO, and H i emissivities. Typical lcr uncertainties in the [C n], 
Hi, and 12 CO emissivities are 0.02, 0.6, and 0.1 Kkms _1 kpc" 1 , 
respectively. We see that both [C n] and 12 CO are mostly con- 
centrated in the 4-10 kpc range, peaking at 4.25 kpc. The dis- 
tribution of the H i emissivities extends over a larger range of 
Galactocentric distances, with a peak at 5 kpc. 



Table 1. Definition of Mask Regions 
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y-. Emission is detected in individual spaxels. 

x: Emission is not detected. 

s/x: Emission is either detected or not detected. 
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Fig. 8. Position velocity maps of the different types of regions defined in Section I47T1 at b — 0°. Mask (grey) represents velocity 
components with only H i detected, Mask 1 (white) are components with only H i and CO detected, Mask 2 (blue) components with 
only Hi and [Cn], Mask 3 (black) components with Hi, [Cn], and CO, and Mask 4 (red) components with Hi, [Cn], 12 CO, and 
13 CO. 



4. The origin of the [Cn] emission in the Galaxy 

4. ) . Velocity Components 

In the following, we associate the [Cn] emission with other 
tracers in order to separate the contribution from different ISM 
phases to the observed emission. We define 5 mask regions 
that represent different stages of the ISM evolution depending 
whether Hi, [Cn], 12 CO, and 13 CO are detected in an individual 
pixel in the position-velocity map (spaxel). We define MaskO 
as those spaxels in which Hi emission is detected but [Cn], 
CO, and 13 CO not. These regions likely represent low-volume 
density atomic gas (see Section l5TI ). Mask 1 includes spaxels in 
which H i and 12 CO (or 13 CO) are detected but [C n] is not. These 
spaxels trace dense molecular gas that is too cold for [C n] to be 
detectable at the sensitivity limit of the observations. Mask 2 are 
spaxels where H i and [C n] are detected but 12 CO is not. These 
regions are likely larger volume density atomic clouds that have 
insufficient FUV-photon shielding to allow the build-up of CO 
and therefore in which most of the gas-phase carbon is in the 
form of C + . Mask 3 includes spaxels in which Hi, [Cn], and 
12 CO are detected but 13 CO is not. These regions represent more 
shielded regions where CO is forming but with inadequate col- 



umn density to have detectable 13 CO emission at the sensitivity 
of our observations. Finally, Mask 4 includes spaxels in which 
Hi, [Cn], 12 CO, and 13 CO are all detected. These regions likely 
represent PDRs that have large volume and column densities 
and are warm enough to produce significant [C n] emission. In 
Table[T]we present a summary of the Mask region definitions. 

In Figure[8]we present the distribution of the different Masks 
in the position-velocity map of the Galaxy for -60° < / <60° 
and b - 0°. We present similar position-velocity maps for 
b = ±0.5° and ±1.0° in Appendix IB] (Figure lB3ll. The dense 
PDRs (Mask 4) are closely associated with the spiral arms. These 
regions are surrounded with spaxels in Mask 3 which, in turn, 
are surrounded with spaxels in Mask 2. We suggest that this spa- 
tial arrangement indicates that Mask 3 and Mask 2 represent the 
envelopes of the active star-forming regions traced by Mask 4. 
There is a significant fraction of spaxels in Maskl, particu- 
larly in the molecular ring and local arms. Mask represents the 
largest fraction of spaxels. 

Figure [9] presents the intensity distribution of Hi, [Cn], 
12 CO, and 13 CO for the different Mask regions defined above 
for b — 0° . Although by definition Masks and 1 have no [C n] 
emission in a given spaxel, its emission can still be seen when the 
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Galactocentric Distance [kpc] 

Fig. 9. Galactocentric distribution of the different types of re- 
gions defined in Section 14.11 The plots include data from all 
LOSs at b = 0°, b = ±0.5°, and b = +1.0°. 



data is averaged over the different rings. The same result happens 
with 13 CO in Mask 3. The detection of emission after averaging 
the data in azimuth shows that our classification into Masks de- 
pends on the sensitivity of our observations. We will study the 
[C n] emission in Mask and 1 and the 13 CO emission in Mask 3 



separately as they represent a different population of clouds that 
are weak in [Cn] and 13 CO, respectively. 



4.2. Atomic gas 

We first assumed that all the observed [C n] emission arises from 
a purely atomic gas at a kinetic temperature 7\j n = 100 K. For 
optically thin emission, the [Cn] intensity (in units of Kkms -1 ) 



-2 



is related to the C + column density, Nc+ ( cm 
den sity of the collisional p artner, n (e~, H, or H2 
e.g. lGoldsmithet~aTll2012l) . 



cm 



and volume 
~ 3 ), as (see 



^[CII] = Nc* 



3.05 xl0 15 (l +0.5(l + i±U 9U1/T * 
R u in 



(1) 



where A u \ - 2.3 x 10 s _1 is the Einstein spontaneous decay 
rate and R u i is the collisional de-excitation rate coefficient at a 
kinetic temperature 7\j n . We used a value of R u i for collisions 
with H at 7\j n = lOOKof 8.1x10' 10 s' 1 cm' 3 (lLaunav & Rouefj 
Il977h . In the optically thin limit we can estimate the atomic gas 
column density, N(H), from the Hi 21cm observations using, 
N(H) = 1.82 x 10 18 /(HI)cirT 2 , with /(HI) in units of Kkms" 1 . 
Observations of H 1 absorption toward continuum sources, how- 
ever, suggest that the opacity of H 1 increases toward the inner 
Galaxy, where the assumption of optically thin emission might 
not apply. iKolpak et all J2002h derived an average value, assum- 
ing Tkin = 50 K, of t =s 1 for R ga \ < 8.5 kpc. For our assumed 
Tkin - 100 K, this opacity corresponds to r =* 0.5. We used 
this opacity to correct the H 1 column densities derived from the 
21 cm line that is associated with [Cn] emission, resulting in a 
30% increase in the derived column density. 

The column density of ionized carbon can be estimated 
from N(Hi) assuming an appropriate [C]/[H] fractional abun- 
dance, provided that all gas-phase carbon is in the form of 
C + . There is evidence that the gas-phase abundance of met- 
als in the Milky Way decreases with Galactocentric distance. 
iRolleston et alj d2000t) found that the fractional abundance dis- 
tribution of several light elements (C, O, Mg & Si) can be rep- 
resented by a linear function with slope -0.07 + 0.01 dex kpc" 1 
over 6 kpc < R sa i < 18 kpc . We use this relationship to convert 
N(H) to N(C + ) for any given galactocentic distance, assuming 
that the slope of the relative abundance gradient does not change 
significantly for R gli i < 6 kpc. We assume a [C]/[H] fractional 
abundance of 1 .4 x°10" 4 for J? ga] = 8.5 kpc (ICardelliet alj|1996l) . 
The [C]/[H] abundance as a function of Galactocentric distance 
is thus given by 



[C]/[H] = 5.5 x 10 4 10 



-0.07« gal 



(2) 



The [C]/[H] abundance drops by a factor of 2.6 from R glL \ - 4 
kpc to 10 kpc. 

With the estimated value of N(C + ) we used Equation ([TJ to 
solve for the volume density of H 1 that is required to match 
the observed [C 11] emission, under the assumption that all ob- 
served [C n] emission comes from atomic gas. In Figure [TOk we 
present the distribution of hydrogen volume densities as a func- 
tion of Galactocentric distance for the different Mask regions for 
4 kpc < R ga \ < 1 1 kpc, which is the range of R gli \ in which most of 
the [C 11] emission is detected. The H volume density increases 
toward the Galactic center with a range between lOOcrrT 3 to 
10 3 cm -3 for Masks 2, 3 and 4. For Masks and 1 the vol- 
ume densities are more moderate, ranging from 20 to 80cirT 3 . 
In Figure [TOb . we present the corresponding thermal pressures 
(p/k - nTkin) as a function of Galactocentric distance for the 



Pineda, Langer, Velusamy, and Goldsmith: A Herschel [C n] Galactic Plane Survey 



3.5 



3 - 



2.5 - 



E 
o 



o 



MaskO — i- 
Mask 1 — *- 
Mask 2 ■■■■»■ 
Mask 3 Ei- 
Mask 4 ■■•■■■■ 



Wolfire et al. 2003; T mi kjn =100 K 



Q.,.Q 



A *..••■. 






.••a .•• 

* 




4 6 8 

Galactocentric Distance [kpc] 



one order of magnitude larger than those observed and predicted 
for the CNM. We therefore conclude that the assumption that 
the [Cn] emission in Masks 2, 3, and 4 originates entirely from 
atomic gas is incompatible with our observations. As we will 
discuss below, part of the observed [C n] emission in Masks 2, 3, 
and 4 also emerges from molecular, and to a lesser extent, from 
ionized gas. 

Motivated by the agreement between the pressures derived in 
Masks and 1 and those predicted in Equation (0, we used this 
expression to estimate, assuming a constant 7\j n = 100 K, the 
radial H volume density distribution the atomic gas associated 
with Masks 2, 3, and 4. We then used Equation ([TJ to estimate 
the contribution from CNM gas to the observed [C n] emission 
in these Mask regions. Note that for a given thermal pressure, the 
intensity of [C n] associated with H i has a only weak dependence 
on the assumed kinetic temperature, varying by about +20% in 
the range between 60 and 200 K. 

In Figure QT] we show the resulting contribution from atomic 
gas to the [C n] intensity as a function of the Galactocentric dis- 
tance. We estimate that the contribution from [Cn] emerging 
from the atomic gas component represents 21% of the observed 
intensity. 
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Fig. 10. (Upper panel) Neutral hydrogen volume density as a 
function of Galactocentric distance for the Mask regions defined 
in Section |4~T1 assuming that all [Cn] is associated with atomic 
gas at 7\i n = 100 K. (Lower panel) Corresponding thermal pres- 
sure as a function of Galactocentric distance. 



different Mask region s. We also show the thermal pressure dis- 
tribution suggested bv lWolfire et al.l (120031) . 



p/k= 1.4xl0 4 exp 



-R 



gal 



5.5 kpc 



Kcm 



(3) 



corresponding to the predicted thermal pressure of the CNM gas. 
Masks and 1 ha ve thermal pressures that are consistent with 
those suggested bv lWolfire et al.1 (12003). Thus, the densities and 
thermal pressures estimated in Masks and 1 are consistent with 
their [C n] emission originating from the CNM. Masks 2, 3, and 
4 have volume densities and thermal pressures that are about 



4.3. Ionized gas 

It has been suggested that a significant fraction of the [Cn] 
emission in the Galaxy observed by COBE emerges from the 
Extended Low Density Warm Ionized M edium (ELDWIM; 
iPetuchowski & Bennett! 1 19931: lHeileslll994T) . We estimated the 
contribution to the observed [C n] emission from an ionized gas 
using the Galacti c distribution of the el ectron density used by 
the NE2001 code (ICordes & Laziol2002l) . The NE2001 code as- 
sumes a Galactic distribution of the electron density, including a 
thin and thick disk, spiral arms (four), and Galactic center com- 
ponents. This electron density distribution is further constrained 
using observations of the dispersion measure of pulsars. 

We used the NE2001 model to extract the electron density 
as a function of heliocentric distance for all GOTC+ lines-of- 
sight in the Galactic plane. For each step in heliocentric distance 
we also calculate the column density of electrons and, assuming 
the [C]/[H] gradient (Equation [2j, the column density of C + as- 
sociated with the ionized gas. We estimate the [C n] intensity as 
a function of heliocentric distance for each line-of-sight using 
Equation (Q]), assuming a kinetic temperature of 10 4 K, and ex- 
citation by collisions with electrons (R„i - 4.8 x 10~ 8 s _1 cm 4 ; 
IWilson & Bellll2002h . The term involving the kinetic tempera- 
ture in Equation ([TJ becomes unity for such a large temperature, 
therefore the intensity of [C n] associated with ionized gas is in- 
sensitive to the kinetic temperature. Note that the NE2001 model 
is smooth and can only reproduce the electron density distribu- 
tion over large scales. Finally, we calculated the distribution of 
the [C n] emissivity as a function of Galactocentric distance tak- 
ing the same considerations as described in Appendix A. The 
resulting radial distribution is shown in Figure QT] We estimated 
that the contribution from [C n] arising from ionized gas repre- 
sen ts only 4% of the total observed intensity. 

IVelusamv et al.l d2013l) observed an excess of [C n] emission 
at the tangent velocities towards the Scrutum-Crux spiral arms 
tangent which was interpreted as coming from compressed ion- 
ized gas. They suggested that this compression of ionized gas 
is the result of the passing of a spiral arm density wave. Note, 
however, that this excess can only be detected toward spiral arm 
tangencies due to a projection effect by which the path length 
toward a spiral arm is maximized. We therefore do not expect a 
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significant contribution from this excess to the overall observed 
[Cn] intensity. 

Note that our estimation of the contribution from ionized gas 
to the observed [C n] emission corresponds only to that for lines- 
of-sight in the Galactic plane with b=0°. T he ELDWIM gas is 
known to have a scale height of about 1 kpc dKulkarni & Heilesl 
119871: lRevnoldsiri989l) . which is much larger than that of the 
atomic and molecular gas (~ 100 - 300pc). Thus, the con- 
tribution from ionized gas to the total [C n] luminosity of the 
Galaxy should be larger than what is suggested by the emissiv- 
ities at b - 0°. At a distance from the Sun of 4 kpc, a 1 kpc 
scale height corresponds to 14°. Thus, within the a 7° beam of 
COBE, the [C n] emission from the ionized component fills the 
beam, while the e mission from the Galact ic plane (vertical an- 
gular scale ~2 deg; iNakagawa et al.|[l998l) will be significantly 
beam diluted. The [C n] emission from the Galactic plane and 
that from the extended ionized medium should appear compa- 
rable within a COBE beam, as suggested b y several analyses 
dPetuchowski & Bennedll993l: lHeiledll994h . A detailed analy- 
sis of the relationship between the Herschel/HIFI [C n] observa- 
tions and those from COBE/FIRAS will be presented in a future 
paper. 

4.4. Dense photon-dominated regions 

It is possible that a significant fraction of the observed [C n] in- 
tensity in the Galactic plane comes from dense PDRs, as obser- 
vations of Galactic and extragalactic PDRs h ave shown that they 
are the sources of b right [Cn] emission (e.g. lStacev et alJfl99U 
iBoreiko et alj|1988l) . Note that the main-beam temperatures in 
the GOTC+ survey are mostly in the 0.5 -8 K range, with only a 
few spaxels with T m b >8 K (see Fig. 7 in lGoldsmith et all2012f) . 
This range of main beam temperatures is well below that ob- 
served in nearby star-forming regions such as Orion (~100K) 
which are powered by the intense FUV radiation field produced 
by nearby massive stars. Bright PDRs are, however, not numer- 
ous enough to dominate the global [Ctrl emissio n in the Milky 
Way. In Section [6] (see also iPineda et al.1 1201 0a) we found that 
most of the [Cn] emission in the Galactic plane arises from 
regions which are exposed to moderate-strength FUV fields, 
in agreement with m odels of the Galaxy based on COBE data 
dCubick et alj|2008l) . We considered [C n] velocity components 
associated with 13 CO (Mask 4) as likely to be associated with 
dense and warm PDRs. Such regions have large enough column 
densities, and presumably sufficient volume densities, to pro- 
duce detectable 13 CO emission while they are also warm enough 
to produce detectable [C n] emission. We determined the con- 
tribution from dense PDRs to the observed [C n] intensity in 
Mask 4 by subtracting the [C n] emission associated with atomic 
gas, as determined in Section 14.21 The resulting radial distribu- 
tion of the [C n] emissivity associated with dense PDRs is shown 
as a red-dashed line in Figure Q~T] The [C n] emission emerging 
from dense PDRs represents 47% of the total emission observed 
in the GOT C+ survey. 

4.5. [Cn] Excess 

The combined contribution to the [Cn] emission from atomic 
gas, ionized gas, and dense PDRs, as estimated above, can ac- 
count for only 72% of the observed emission. It has been sug- 
gested the remaining [C n] emission is produced in the envelopes 
of dense molecular clouds where the H 1/H2 transition is largely 
complete, but the column densities are low enough that the 



5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 
0.5 





Total [CM] Observed 
[CM] from dense PDRs 
[CM] from Atomic Gas 
[CM] from Ionized Gas (NE2001) 








5 10 

Galactocentric Distance [kpc] 



20 



Fig. 11. Azimuthally-averaged radial distribution of the ob- 
served [C n] emissivity in Masks 2, 3, and 4 together with the 
different contributions from the ISM components discussed in 
Section [4] Typical uncertainties are 0.02 K km s" 1 kpc 1 for the 
total [C 11] emission and that contributed from dense PDRs, and 
0.001 Kkms 1 kpc _1 for the contribution from Hi gas. For the 
contribution from ionized gas, we assume 30% uncertainties or 
0.07 K km s~' kpc~'. 



C + /C °/CO transition is not (e.g lMadden et al.ll997tlLanger etall 
120101) . Since this H2 gas component is not traced by CO, but by 
[C 11] and to a lesser extent by [C 1], we refer to it as "CO-dark 
H2 gas". In the following section we derive the distribution of 
the CO-dark H2 gas in the plane of the Milky Way, and compare 
it to that traced by 12 CO and 13 CO. 

5. The distribution of the ISM components in the 
plane of the Milky Way 

5. 1 . The Warm and Cold Neutral Medium 

As mentioned earlier, atomic gas is predicted to be in two 
phases in roughly thermal pressure equilibrium, the Cold Neutral 
Medium (CNM; « H =* 50 cm -3 ; 7/ kin ^80 K) and the Warm 
Neutral Medium (WNM; n H =* 0.5ctrT 3 ; 7\ in ^8000 K). 
Because atomic hydrogen is one of the main collisional part- 
ners of C + , the [Cn] line is a tracer of atomic gas that has the 
advantage over the 21 cm line (observed in emission) that it is 
sensitive to both the volume density and the kinetic temperature 
of the gas (Equation[TJ. Thus, the [C 11] line should in principle 
trace both the WNM and CNM atomic gas components. In prac- 
tice, however, because of the volume density contrast between 
the CNM and WNM, the [C n] line is expected to be significantly 
weaker for WNM ga s compared to the emission from CNM gas 
dWolfireet all 120101) . For the Hi peak intensity of 160K (in a 
0.8 km s -1 channel width) at the 5 kpc ring, and considering the 
abundance and thermal pressure gradients (Equation |2]and [3] re- 
spectively), we can use Equation ©, assuming 100K for the 
CNM and 8000 K for the WNM, to estimate a [C n] intensity as- 
sociated with H 1 of 0.3K in the case that all gas is CNM, and of 
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Fig. 12. Radial distribution of the azimuthally averaged total H i 
column density in the Galactic plane. We also include the es- 
timated radial distribution of the cold neutral medium (CNM) 
and warm neutral medium (WNM) gas. Typical uncertainties of 
the Hi column densities are 10 18 cm~ 2 kpc _I . The vertical line 
shows the Solar radius. At 9 kpc there is a discontinuity in the 
total and WNM JV(Hi) distribution which is attributed to the 
abrupt change between the number of rings sampled in the inner 
and outer Galaxy. We therefore also show lines that smoothly 
connect the distributions in the inner and outer Galaxy. 



0.04K in case all gas is WNM. Considering our 3<x sensitivity 
limit per channel of 0.3 K we find that weak CNM gas can be 
detected, while WNM gas would be too weak to be detected in 
individual spaxels. 

When azimuthally averaging (Section 14. 2t , however, [C n] 
emission associated with CNM can produce a significant con- 
tribution to the observed [C n] emissivity. For the R gi i - 
5kpc peak of 1.4 x 10 21 cm _2 kpc _1 , we find a contribu- 
tion of 1.1 K km s 1 kpc 1 for all Hi gas in CNM, and 
O.CWKkms'kpc 1 for all Hi gas in WNM. As the 3cr sen- 
sitivity limit for the azimuthal average is 0.06 K km s 1 kpc -1 
(AppendixlAl. we see that WNM is always below the sensitivity 
limit even for the azimuthal average. This suggests that only the 
CNM can contribute to the observed [C n] emission. We there- 
fore assume that spaxels that have both H i and [C n] emission 
(Masks 2, 3, and 4) are associated with CNM gas. In our survey, 
there are also spaxels where H i and 12 CO emission are detected 
but [Cn] emission is not (Mask 1). These regions are likely too 
cold and/or have too low volume densities for [C n] to be de- 
tected, but can still be associated with cold H i gas. We therefore 
also consider H i associated with Mask 1 as being in the CNM 
phase. Finally, we consider spaxels with H i but neither [C n] nor 
12 CO emission (MaskO) to be associated with WNM gas. In the 
following we use this criteria to separate the CNM and WNM 
components from the observed H i emission and study their dis- 
tribution in the plane of the Milky Way. 

It is possible that some spaxels have [C n] emission that is 
below the sensitivity limit for individual spaxels but can still be 
detected in the azimuthal average. In Figure [9] we see that [C n] 
emission (~0.1Kkms~' kpc -1 ) in the inner galaxy can be de- 
tected in the region where only H i is detected in the pointed ob- 
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Fig. 13. Radial distribution of the mass fraction of the atomic gas 
in the cold neutral medium (CNM). The hori zontal line shows 
the loc al mass fraction of the CNM derived bv lHeiles & Trolandl 
d2003l) . The vertical line denotes the solar radius. The typical 
uncertainty in the ratio is ~0.002. 



servations (Mask 0). The average H i column for R gii \ = 3-8 kpc 
(~ 1 x 10 20 cm _2 kpc _1 ), would imply a [Cn] emissivity of 
about 0.17Kkms _1 kpc -1 in the case that all gas is CNM and 
0.006 Kkms" 1 kpc _1 in the case that all gas is WNM. Thus, the 
emission detected in Mask can only be associated with CNM 
gas. Note that for the inner galaxy, when we associate H i with 
[C n] and/or 12 CO emission in individual spaxels, we are assum- 
ing that the H i and [C n] and/or 12 CO emission arises from both 
the near and far side of the Galaxy. We could in principle over- 
estimate the H i column density associated with the CNM in the 
case that the [C n] and/or 12 CO emission arises from either the 
near or far side of the galaxy, but H i from both sides. We con- 
sider this situation unlikely, however, as we expect that the phys- 
ical conditions within a galactocentric ring to be relatively uni- 
form. The maximum effect, in the case that all LOSs in the inner 
galaxy are affected, is a reduction of the azimuthally averaged 
CNM column density by a factor of ~2 (and a corresponding 
increase of the WNM column density by the same factor). 

In Figure Q~2] we show the radial distribution of the total 
H i column density together with that of the contributions from 
CNM and WNM gas. The column densities are estimated from 
the H i emissivity distributions which are in turn calculated con- 
sidering all spaxels with Hi emission as well as those associ- 
ated with CNM gas (Masks 1, 2, 3, and 4) and those associ- 
ated with WNM gas (MaskO). The column density of CNM gas 
was corrected for opacity effects as described in Section 14.21 
The value of iV(Hi) at R gai = 8.5 kpc (6.7 x 10 20 cm- 2 kpc _1 ; 
without the opacity correction) is in good agreement with the lo- 
cal value of the H i surface density of 6.2 x 10 20 cirT 2 found by 
iDickev & Lockmanld 19901) . As we can see, the CNM dominates 
the H i column density in the inner Galaxy, while WNM domi- 
nates in the outer Galaxy. The column density of the WNM gas 
in the inner Galaxy is similar to what we estimated for the H i 
column density required to reproduce the [C n] emission result- 
ing from the azimuthal average in Mask 0. This result indicates 
that even a larger fraction of the atomic gas in the inner Galaxy 
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might be in the form of CNM. CNM clouds may still exist in the 
outer Galaxy, as molecular clouds and star formation are present, 
and the reduction of its contribution in the outer Galaxy might 
be a result of a reduced filling factor of the CNM gas. 

There is a discontinuity in the total H i column density distri- 
bution at around 9 kpc, which may be an artifact of the assumed 
flat rotation curve with purely circular motions. An alternative 
explanation is shown in Figure lA72l in Appendix|A] which shows 
the number of ring samples as a function of Galactocentric dis- 
tance. Because a LOS samples a ring in the inner Galaxy twice 
but a ring in the outer Galaxy only once, there is a discontinuity 
in the number of samples between the inner and outer Galaxy. 
This discontinuity will therefore affect the azimuthally averaged 
emissivity distributions. The effect on the [Cn] and 12 CO emis- 
sivity distributions is small because it coincides with a steep de- 
crease in the value of the sum of the emissivities for each sam- 
ple (see numerator in Equation lA.8b . In the case of H i, however, 
the sum of the emissivities for each sample is smooth at around 
9 kpc, and therefore dividing by the number of samples results 
in a noticeable discontinuity. In Figure [12] we also show the Hi 
column density distribution where data points with R„ d \ < 9 kpc 
and R gSL \ > 1 1 kpc are smoothly connected. 

In Figure Q~3] we present the mass fraction of the CNM 
gas. In the inner Galaxy, up to 90% of the atomic gas is in 
the form of CNM, with even a larger percentage being pos- 
sible based on the discussion above. For larger Galactocentric 
distances, this fraction decreases to only 0.1-0.2, showing lit- 
tle variation with Galactocentric dista nce. These re s ults ar e in 
good agreement with those found by iDickev et al.1 (l2009h us- 
ing Hi emission/absorption observations in the outer Galaxy. 
At R gd \ - 8.5 kpc, we find that the fraction is in good agree- 
ment w ith the local fraction of 0.4 derived bv lHeiles & Trolandl 
(120031) . We derive an average CNM mass fraction for the entire 
Galactic disk of 0.43. 

Our results suggest that the CNM is the dominant atomic gas 
phase in the inner Galaxy while the WNM dominates the outer 
Galaxy. The increase of the fraction of CNM gas towards the 
inner Galaxy coincides with that of the 12 CO and [C n] emissiv- 
ities. This suggests that the CNM plays a important role as the 
precursor of molecular clouds that will eventually form stars. 

5.2. Molecular Hydrogen 

5.2.1. CO-dark H 2 gas 

To derive the column density radial distribution of the CO-dark 
H2 gas component, we first estimated the C + column density as- 
sociated with the [C 11] excess emission derived in Section |4] 
The C + column density, which can be derived by inverting 
Equation (Q]), depends on the observed [Cn] intensity, the H2 
volume density, «h, , and the kinetic temperature of the H2 gas, 
7^. Without an independent estimate of the physical conditions 
in the CO-dark H2 layer, we need to make an assumption about 



the value of «h, and 7^. In Section 14721 we derived the con- 
tribution from atomic gas (CNM) to the observed [C 11] inten- 
sity by estimating the H 1 volume density, rim, from the radial 
thermal pressure distribution (Equation |3), assuming a kinetic 
temperature of the Hi gas, T™, of 100 K. After the H/H 2 tran- 
sition takes place, and as the column density increases, self- 
gravity starts to become important, resulting in a gradual in- 
crease of the volume density of the gas. The increased shield- 
ing of FUV photons results in a decreased gas heating, while 
line cooling is largely unaffected due to the larger volume densi- 
ties, resulting in a reduction of the kinetic temperature. Overall, 
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Fig. 14. Radial distribution of the azimuthally averaged excess 
[C n] emissivity as discussed in Section l4~5l ( upper panel). Radial 
distribution of the C + column density derived from the excess 
[C 11] emissivity (middle panel). Radial distribution of the H2 
column densities derived from N(C + ) associated with the excess 
[Cn] emissivity (lower pannel). 



however, the thermal pressure of the gas increases with column 
density: studies of local clou ds, in CO absorption in the mil- 
limeter (lLiszt&Lucaslll998l) and in the UV (Goldsmith 2013 
submitted), tracing column densities up to that of the C + /C°/CO 
transition layer, typically suggest thermal pressures in the range 
between 10 2 to 10 4 Kcm~ 3 , while for well shielded regions, 
traced by 13 CO, excitation studies su ggest thermal pressu res in 
the range between 10 4 to 10 5 Kcm" 3 (ISanderset al.lll993l) . It is 
therefore likely that a thermal pressure gradient is present in the 
region between the H/H2 and C + /C°/CO transition layers, where 
the emission from CO-dark H2 gas originates. To estimate the 
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Fig. 15. CO-dark H2 gas column density as a function of 
Galactocentric distance for different Hi and H2 temperature 
pairs (T™;?^). The corresponding volume densities were de- 
rived from Equation (|3] see Section l5.2.U . 



column density of CO-dark H2 gas, we assume that the thermal 
pressure of the CO-dark H2 gas is the geometric mean between 
that of the Hi gas (Equation [3} and the maximum observed at 
the C + /C°/CO transition layer. At the solar neighborhood, the 
minimum pressure given by Equation [3] is 3000 K crrT 3 while 
the maximum pressure at the C + /C°/CO transition layer given 
by Goldsmith (2013) is 7400 KcrrT 3 . The resulting geometric 
mean is 4700 KcrrT 3 , which is a factor of 1.6 larger than the 
thermal pressure given by Equation ((3}- We therefore multiply 
Equation 0} by this factor to determine the mean thermal pres- 
sure of the CO-dark H2 gas as a function of Galactocentric dis- 
tance. We assume a consta nt kinetic temp e rature of T. 2 =70K. 



following the results from iRachford et aL (120021). and a corre- 
sponding R u , = 4.3 x 10' 1U s' 1 cm 3 dFlower & Launavlll977l) 
for collisions with H2. Using this kinetic temperature we derive 
the H2 volume density of the CO-dark H2 gas from Equation © 
as discussed above. We will discuss below the effects of choos- 
ing different values of 7*^ and «(H2) in the derivation of the 
column density of the CO-dark H2 gas. We finally converted the 
radial distribution of the C + column density to that of the H2 col- 
umn density, assuming that all gas-phase carbon is in the form 
of C + , by applying a [C]/[H2] abundance gradient. The [C]/[H2] 
abundance ratio is twice the [C]/[H] ratio given in Equation @. 
In Figure [14] we show the azimuthally-averaged radial dis- 
tribution of the [C 11] excess emissivity (top panel), of the corre- 
sponding C + column density (middle panel), and of the CO-dark 
H2 column density (lower panel). The steep distribution of the 
[Cn] excess emissivity, with a peak at 4.2 kpc, becomes flatter 
when converted to N(C + ) as a result of applying the «(H2) gradi- 
ent derived from Equation (01 and a constant T^ equal to 70 K. 
When applying the [C]/[H2] abundance gradient, the radial dis- 
tribution of CO-dark H2 becomes nearly constant between 4 and 
1 1 kpc. We will discuss this flat distribution of the CO-dark H2 
column density and relate it to the denser H2 gas traced by CO 
and 13 CO in Section l531 Note that when calculating column den- 



sities we are assuming that the [C 11] emission is extended with 
respect to the beam of our observations. In the case that there is 
substructure within the beam, our estimate of the column den- 
sity will be a lower limit, as it corresponds to that averaged over 
HerscheFs 12" beam. 

In Figure[15]we show the resulting CO-dark H2 column den- 
sity distribution for different combinations of the kinetic temper- 
ature in the H 1 and H2 layers. The volume densities are obtained 
from Equation (fJJ, considering the factor 1.6 correction for H2 
gas discussed above. We assume that the temperature of the H2 
layer is equal to or lower than the temperature of the H 1 layer. 
For each kinetic temperature, the corre sponding value of R„i was 
interpolated from data calcul ated by Flower & Launavl (119771) 
for collisions with H2 and by lLaunav & Rouefil (1 19771) for col- 
lisions with H. As we can see, for a wide range of temperature 
combinations, the derived column density deviates by about 30% 
from our adopted model with 7"™ = 100K and 7^ = 70 K. We 
therefore consider that our CO-dark H2 column density deter- 
mination has a 30% uncertainty. 

5.2.2. H 2 gas traced by I2 CO and I3 CO 

A commonly employed method to estimate the H2 mass of 
molecular clouds is the use of the 13 CO integrated intensity, 
7( 13 CO), to derive the 13 CO column density, JV( I3 CO), assuming 
that this line is optically thin and that local thermodynamic equi- 
librium (LTE) applies. JV( 13 CO) can be then converted to N(CO), 
assuming an appropriate isotopic ratio, and N(CO) to iV(H2), as- 
suming a CO abundance relative to H2. It has been observed in 
large-scale 12 CO and 13 CO surveys, however, that the spatial ex- 
tent of the I3 CO emission is often noticeably smaller than, and 
is contained within, the observed extent of the 12 CO emission 
(iGoldsmith et al.l 120081; iRoman-Duval et al.ll2010l) . This spatial 
discrepancy is mostly due to the limited sensitivity of the obser- 
vations that is often inadequate to detect the weak 13 CO emission 
emerging from the envelopes of molecular clouds. The region 
where only 12 CO is detected but 13 CO is not can account for 
~30% of the total mass o f a molecular cloud dGoldsmith et al.l 
120081: iPineda et al1l2010al) . The discrepancy between the spatial 
extent of these two lines is also seen in the position velocity 
maps of the observed 12 CO and 13 CO lines obtained as part of 
the GOT C+ survey. 

To estimate the azimuthally averaged H2 column density 
traced by CO and 13 CO for each ring, we considered two cases: 
spaxels where both 12 CO and 13 CO are detected (Mask C02) and 
spaxels where 12 CO is detected but 13 CO is not (Mask COl). In 
FigureQ~6] we show the radial distribution of the azimuthally av- 
eraged I2 CO and 13 CO emissivity for spaxels in Masks COl and 
C02. Although by definition spaxels in Mask COl have 13 CO in- 
tensities below the sensitivity limit, we could nevertheless detect 
I3 CO emission for spaxels in this mask region after averaging 
them in azimuth as seen in the figure. As 13 CO can be considered 
optically thin, the integrated intensity of this line is proportional 
to the 13 CO column density which in turn should be proportional 
to the H2 column density. Thus, spaxels in Masks C02 and COl 
correspond to regions with larger and smaller column densities, 
respectively, with the threshold between these two column den- 
sity regimes given by the sensitivity of the 13 CO observations. 

In the following we describe our procedure used to estimate 
H2 column densities from Mask COl and MaskC02. The pro- 
ce dure is very similar to that used in the Tauru s molec ular cloud 
by IPineda et all d2010al) and IGoldsmith et"ai] d2008l) . We refer 
to these papers for further details on the column density de- 
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termination. We derived the H2 column density using the 12 CO 
and 13 CO emissivities calculated from spaxels where the 12 CO 
and/or 13 CO are detected. We later expressed the derived H2 col- 
umn density as an average that considers all sampled spaxels 
within a given Galactocentric ring, not only where 12 CO and 
13 CO emission is detected, to facilitate the comparison with the 
column density of CO-darkH2 gas (Section[5]2j}. 

For MaskC02, we estimated N( 13 CO) from 7( 13 CO) under 
the assumption that this line is (almost) optically thin and that 
LTE applies. Apart from the dependence on 7( 13 CO), JV( 13 CO) 
depends on the excitation temperature, T ex , and opacity correc- 
tions that might apply. The excitation temperature and the optical 
depth can be derived from the peak line intensity of the 12 CO and 
13 CO lines, respectively. To derive the line peak intensities we 
assumed that the average 12 CO or I3 CO emissivities are equiv- 
alent to the average of the integrated intensity of a number of 
Gaussian lines with equal FWHM line width of 3 km s -1 , which 
is a typical value for the observed 12 CO and 13 CO lines in our 
survey (Langer et al. 2013 in preparation). Note that the calcu- 
lation of r ex also assumes that the 12 CO line is optically thick 
(r >> 1). Accounting for opacity in the derivation of JV( 13 CO) 
typically results in only a 10% correction. 

We converted JV( I3 CO) to JV(CO) us ing the slope of th e 
12 C/ 13 C isotope ratio gradient derived by ISavage et al.l (120021) . 
while fixing the 12 C/ 13 C isotope ratio to be 65 at R ga \ - 8.5 kpc. 
The 12 C/ 13 C isotope ratio as a function of Galactocentric dis- 
tance is thus given by 



13 C 



71^1- 
v kpc/ 



(4) 



where R g!L i is the Galactocentric distance. For this gradient, the 
isotopic ratio varies from ~43 to ~80 for Galactocentric dis- 
tance between 4 kpc and 12 kpc. We converted N(CO) to N(H2) 
using a [CO]/[H2] abundance gradient constructed using the 
slope of the [C]/[H] abundance gradient in Equation © fixing 
[CO]/[H 2 ] = l x 10 - 4 atR m] =8.5 k pc. 

As noted by e.g. lTielensl(l2005l) . the [CO]/[H 2 ] relative abun- 
dance used here, and widely used in the literature, is a fac- 
tor of 3 lower than the gas-phase carbon abundance observed 
in the diffuse ISM, used in Section 15.2.11 This reduced CO 
fractional abundance is inconsistent with chemical models that 
predict that, in well shielded re gions, all the gas-phase car- 
bon should be in the form of CO (ISolomon & Klempererjll972t 
iHerbst & Klempererlll973l) . This discrepancy is not a result of 
CO molecules frozen onto dust grains, as it persist even af- 
ter a correction to account for CO an d CO2 ices is applied 
dWhittet et al.l2010HPineda et al.l2010al) . The CO/H 2 abundance 
is often derived by comparing observations of CO and 13 CO with 
visu al extinction maps towards well studied nearby dark cloud s 
(e.g. lDickman|[l978t iFrerking et al.lll982l iPineda et al.ll2010al) . 
When the contribution from H 1 gas is subtracted, the visual ex- 
tinction traces the total column of H2 along the line of sight. 
However, a certain column of H2 is required to shield the FUV 
photons before CO can form. Therefore, for a given path length, 
only a fraction of it, its central part, will have all gas-phase car- 
bon in the form of CO, while the remainder will be what we 
call the CO-dark H 2 gas. Dividing the column density of CO 
with that of H2 will result in a fractional abundance that is lower 
than what would result if CO is the dominant form of carbon 
over the entire path length. This effect is more prominent in re- 
gions with moderate column densities, which represent the bulk 
of the mass in giant molecular clouds, where the CO-dark H2 
and CO-traced layers are comparable. Following this argument, 
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Fig. 16. Radial distribution of the azimuthally averaged 12 CO 
and 13 CO emissivities observed in the Galactic plane for the 
case where 12 CO is detected in individual spaxels but 13 CO is 
not (Mask COl ; upper panel), and for the case where both 12 CO 
and 13 CO are detected in individual spaxels (Mask C02; lower 
panel). Typical uncertainties of the I2 CO and I3 CO emissivities 
are 0. 1 and 0.02 K km s kpc -1 , respectively. 



the local values, [C]/[H 2 ]=2.8 x 10~ 4 and [CO]/[H 2 ] = l x 10~ 4 , 
would imply a CO-dark H 2 gas fraction of 0.6, which is con- 
sistent w ith what is shown in Figure [19] and that independently 
found bv lParadis et al.l d2012l) studying the CO-dark H2 fraction 
in the Solar neighborhood using extinction mapping. Note that 
by using [CO]/[H2]=l x 10~ 4 in our calculations we are implic- 
itly accounting for the contribution from CO-dark H2 gas in the 
derivation of the H2 column density from CO and 13 CO The 
CO-dark H2 gas present in regions traced by CO and 13 CO will 
be studied in detail by Langer et al. (2013 in preparation). 

As noted above, spaxels in Mask COl correspond to low H2 
column density gas. For molecular clouds it is likely that this 
gas is also characterized by low IT volume densities. Assuming 
that the detected '"CO emission in Mask COl is nonetheless op- 
tically thick we estimated the excitation temperature as we did 
for Mask C02. Typical values were found to be in the 4K to 
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6K range for 4kpc < R sa \ < 12kpc. As seen in Taurus, pixels 
in Mask COl are located in the periphery of the large column 
density molecular cloud traced by I2 CO and 13 CO (MaskC02). 
However, the excitation temperatures estimated for Mask COl 
are lower than those derived for Mask C02 (~ 10-20 K), and 
therefore they are inconsistent with what we would expect for 
this gas being in LTE, because the kinetic temperature is likely 
to be higher in the outer regions of the cloud which are subject 
to increased heating. It is therefore reasona ble to assume that the 
gas in Mask COl i s subthermally excited. Pineda et al.1 (1201 Oal) 
used the RADEX (Ivan der Tak et al.ll2007l) excitation/radiative 
transfer code to predict line intensities from a gas with condi- 
tions similar to those of Mask COl in the envelope of the Taurus 
molecular cloud. RADEX accounts for effects of trapping, which 
are important for the excitation of CO at low densities. The 
model parameters are the kinetic temperature, T^ n , the H2 vol- 
ume density, ^(IT), the CO column density per unit line width, 
N(CO)/6v, and the 12 CO/ 13 CO abundance ratio, R. The observed 
parameters are the excitation temperature and the 12 CO/ 13 CO in- 
tegrated intensity ratio. As mentioned above, the excitation tem- 
perature varies from 6K to 4K for 4 kpc < /? ga i < 12kpc. The 
12 CO/ 13 CO integrated intensity ratio shows little variation from 
its average value of ~24 for the same range in Galactocentric 
distance. We assume that the kinetic temperature is 15 K. If the 
value of R were known, one could find a unique solution for 
n(H2) and N(CO). However, the isotopic ratio is expected to 
vary from that expected for well-shielded regions due to iso- 
topic enhancements produced by chemical and/or photo effects. 
Thus, there is a f amily of solutions th at can reproduce the ob- 
served quantities. iPineda et al.1 (1201 Oal) used nearly 10 6 individ- 
ual spectra in the Taurus molecular cloud to find a set of solu- 
tions that produced a monotonically decreasing R with decreas- 
ing excitation temperature and a smoothly decreasing column 
density with decreasing excitation temperature. The observed 
12 CO/ 13 CO integrated intensity ratio for a given value of T ex 
for all Galactocentric rings is in excellent agreement with the 
same quantities observed in the envelope of the Taur us molecular 
cloud. W e therefore used the solution presented in IPineda et alJ 
(l2010ah to derive the N(CO) column density in Mask COl for 
each Galactocentric ring. For our observed values of T ex and the 
12 CO/ 13 CO integrated intensity ratio, their solution suggests a 
gas with n(H2) - 250 cirT 3 and R =* 30. We use their fit to the 
relationship between CO column density and the observed exci- 
tation temperature given by 
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to derive the average N(CO) column density in Mask COl. 

The [CO]/[H2] abundance ratio is expected to show strong 
differences between UV-exposed and shielded regions as a re- 
sult of the c ompetition between the forma t ion and destructio n 
of CO (e.g Ivan Dishoeck & Blackl Il988t IVisser et al.l 120091) . 
IPineda et alJ d2010al) compared CO column densities derived at 
the envelope of the Taurus molecular cloud with a visual extinc- 
tion, Ay, map of the same region (with the contribution from 
dust associated with atomic gas to the visual extinction sub- 
tracted). They provided a polynomial fit giving Ay at any given 
value of N(CO) (see their Figures 13 and 14). We use this fit to 
derive visual extinctions from N{CO) in Mask COl. Since Ay 
traces the total column density of H2, not only that traced by 
CO, the relationship between N(CO) and Ay should already ac- 
count for the CO-dark H2 gas present in the envelope of the 
Taurus molecular cloud. The visual extinction outside the CO 
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Fig. 17. Radial distribution of the azimuthally averaged H2 col- 
umn density as traced by 12 CO and 13 CO We include the radial 
distribution of A^Hz) from gas associated with spaxels where 
both 12 CO and 13 CO are detected (Mask C02) and where 12 CO 
is detected but 13 CO is not (Mask COl). The total distribution 
(Mask COI+ Mask C02) has been increased by 20% to account 
for 12 CO emission that is not detected in individual spaxels but 
was detected after averaging them in azimuth (see Section l5.2.2l >. 



boundary in Taurus can be estimated by searching for the value 
of Ay at which N(CO) = OctrT 2 . A polynomial fit to the lower 
part of the N(CO)-Ay relationship gives N(CO) = cirT 2 for 
Ay=0.37mag. Note that this correction only applies if the col- 
umn density in Mask COl is derived using spaxels that have 
both [Cn] and 12 CO emission. There are, however, spaxels in 
Mask COl where I2 CO emission is detected but [C 11] is not, and 
the correction in the conversion from N(CO) to Ay does not ap- 
ply to them, because their CO-dark H2 layer is not already traced 
by [C n] . We repeated the column density estimation described 
above but for only spaxels in Mask COl that have both [C n] and 
I2 CO emission as well as for those that have 12 CO emission but 
no [Cn]. After applying the 0.37 mag correction to the column 
density distribution derived from spaxels detected in [C 11] and 
I2 CO and adding it to that derived from spaxels with 12 CO emis- 
sion only, we find that the total column density distribution is 
about 15% lower than that without applying the correction. We 
therefore reduce the derived values of Ay by 15% to account for 
the CO-dark H 2 gas already traced by [C 11] in Mask COl. 

The formation of CO in cloud envelopes is mostly influ- 
enced by dust grains shielding the FUV radiation field, with 
smaller contributions from H2 shielding and CO self-shielding. 
Thus, changes in metallicity will not significantly influence 
the dust column density, but will modify the H2 column den- 
sity, at which the C + /C°/CO transition takes place, because 
a reduction/increase on metallicity is associated with a in- 
crease/reduction of the gas-to-dust ratio, R gl \ = N(H2)/Ay. To 
convert the derived visual extinctions to H2 column densities, we 
assumed that R s a has the same slope, but with the opposite sign, 
of the galactic metallicity gradient (Equation|2), an d we fixed its 
value at i? gal = 8.5 kpc to be 9.4 x 10 20 cnr 2 mag -1 dBohlin et all 
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Il978l) . The resulting gas to dust ratio as a function of galacto- 
centric distance is therefore given by, 



Rod = 2.38 x 10- u 10' 



20 ln 0.07« gal 



cm /mag, 



(6) 



with the galactocentric distance, R gd \, in units of kpc. We applied 
thi s relationship t o conve rt fro m Ay to JV(H2) in Mask CO 1 . 

iPineda et all d2010al) and iGoldsmith et al.1 d2008l) detected 
12 CO and 13 CO emission in the periphery of Taurus where these 
lines were not detected in individual pixels. [Pineda et al1(l2010al) 
estimated that this gas accounts for about 20% of the total mass 
of Taurus. Although we can still detect weak 12 CO emission af- 
ter averaging in azimuth regions where 12 CO is not detected in 
individual spaxels, we could not detect I3 CO, and therefore can- 
not repeat the analysis done in Taurus. We corrected the total 
column density traced by 12 CO and 13 CO by 20% to account for 
this missing CO gas component. 

Figure [17] shows the total azimuthally averaged H2 column 
density as traced by 12 CO and I3 CO For the range between 4 kpc 
and 7 kpc, gas associated with Mask C02 contributes about 70% 
ofthetotalH 2 traced by CO and 13 CO. For7kpc< fl ga i <10kpc, 
the gas associated with Mask C02 contributes about 25% of the 
total CO-traced FF gas. This contribution is consistent with the 
relative contribution from Mask COl and Mask CQ 2 to the mass 
of the n earby Taurus molecular cloud derived by IPineda et al.l 
(I2010ah . 



5.3. The complete distribution of H 2 in the plane of the Milky 
Way 

In Figure [18] we show the radial distribution of the total az- 
imuthally averaged IT column density. The total iV(FF>) is the 
combination of the contribution from H2 gas traced by CO and 
that traced by [Cn] (CO-dark H2 gas). Both contributions are 
also included in the figure. Close to the Galactic center, most 
of the H2 is traced by CO with the CO-dark IT gas making a 
negligible contribution. The CO-traced H2 gas is mostly con- 
centrated in the 4 kpc to 7 kpc range, while the CO-dark FF gas 
is extended over a wider range of Galactocentric distances, be- 
tween 4 and 1 1 kpc. As mentioned in Section 15.2.11 the distri- 
bution of the CO-dark H2 column density is nearly constant be- 
tween 4 and 1 1 kpc. However, the fraction of the total molecular 
gas in this component increases with Galactocentric distance. 
This increase is shown in Figure [19] where we show the CO- 
dark H2 gas fraction as a function of Galactocentric distance. 
The fraction of CO-dark H2 gas increases monotonically with 
Galactocentric distance, rising from 0.2 of the total H2 gas at 
4 kpc to be about 0.8 at 10 kpc. The Galactic metallicity gradi- 
ent is likely accompanied by a decrease of the dust-to-gas ratio 
with Galactocentric distance. The reduced dust extinction results 
in an increased FUV penetration, which in turn results in the 
C + /C°/CO transition layer taking place at larger H2 column den- 
sities. Thus, as suggested by our results, the low-column density 
H2 gas at large Galactocentric distances, with low-metallicities, 
is better traced by [C 11] (and perhaps [C 1]) than by CO. 

For R gd i > 1 1 kpc, we see a reduction of the column density 
of CO-dark H2 gas. In the outer Galaxy we detect fewer [C n] 
clouds per unit area and, as we can see in Figure Q~T] most of 
the detected emission is associated with dense PDRs. Because 
of the reduced carbon abundance in the outer Galaxy, and per- 
haps reduced thermal pressures, the [C 11] intensity associated 
with CO-dark H2, as well as that associated with H 1 and electron 
gas, might be reduced to be below our detection threshold, while 
only regions associated with star-formation, dense PDRs, can 
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Fig. 18. Radial distribution of the H2 column density in the plane 
of the Milky Way. We also show the relative contributions to the 
total H2 column density from gas traced by 12 CO and 13 CO and 
from CO-dark FF gas. 
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Fig. 19. Radial distribution of the fraction of the total FF column 
density comprised by the CO-dark H2 gas. 



produce detectable emission. Due to the reduced star formation 
activity per unit area in the outer Galaxy, the average strength of 
the FUV field is likely also reduced. This reduction of the FUV 
field results in a lower kinetic temperature that reduce the exci- 
tation of [C 11] line. Additionally, the plane of the Milky Way is 
known to warp towards the outer Galaxy, while the CO-dark FF 
distribution is based on [C 11] data at b-Q° . Thus, our sampling of 
the outer Galactic plane might not be as good as it is for the inner 
Galaxy. Therefore we cannot rule out that the distribution of FF 
extends even further in the outer Galaxy than indicated here. 
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Fig. 20. Radial distribution of the CO-to-H2 conversion factor. 
We show the distribution for values estimated considering the 
CO-traced H2 gas as well as that including the CO-dark H2 gas 
component. 

5.4. The CO-to-H 2 conversion factor in the Milky Way 

When 13 CO observations are not available, the mass of molec- 
ular clouds or the molecular content of entire galaxies is of- 
ten estimated from observations of the 12 CO / = 1 — > line 
applying an empirically-derived CO-to-H2 conversion factor 
(Xco = N(H2)/Ico K MuJLco). Modelling applied to y-ray 
observations, which trace the total molecular content along the 
line-of-sight, results in a local value of 1 .74 x 10 2 cm -2 (Kkm 
s" 1 )- ' or M Hl /Lco = 3.7M (K Km s" 1 pc 2 )" 1 dGrenier et al.l 
120051) . Several theoretical studies have focused on the depen- 
dence of Xqo on environmental conditions, showing that Xqo 
is particularly sensitive to metallicity and the stre n gth of the 
FUV radiation field (e.g. iMalonev & Blackl 119881: lSakamotd 
119961) . Observationally, however, it has been found tha t Xqq is 
mostly sensitive to meta llicity (e.g. iRubio et all 1 199 lb I Wilson! 
119951: llsraell 1 1 997J l2000l) while not showin g a significant de- 
pendence on the st rength of the FUV field dPineda et al.ll2009t 
iHughesetalhOlOh . 

We estimated the radial distribution of Xqo in the 
Galactic plane by dividing the H2 column density in a given 
Galactocentric ring by its corresponding 12 CO emissivity. In 
Figure [20] we show the radial distribution of Xco derived from 
the H2 column density traced by CO and 13 CO only, as well 
as that including the CO-dark H2 component. We also show 
the radial distribution of Xqo derived from the relati onsh ip be- 
tween Xqq and metallicity observed bv lWilsonl d 1995b and llsraell 
d2000l) in nearby galaxies. We converted the dependence of Xco 
on metallicity presented by these authors to a radial dependence 
using the slope of the metallicity gradient in Equation d2}, fixing 
X co to be 1.74 x 10 20 cirr 2 (K km s" 1 )- 1 at fl gal = 8.5 kpc. 

Both Xco distributions increase with the Galactocentric dis- 
tance. But when the contribution from CO-dark H2 is included, 
the slope is steeper. This steeper slope is a result of the mono- 
tonic increase of the CO-dark H2 contribution to the total H2 
column density with Galactocentric distance (Figure[T9t. The ra- 
dial distribution of Xqo considering JVfFb.) from 12 C O and 13 CO 
only is consistent with that fitted by IWilsonld 19951) . When CO- 



dark H2 is included in the calculation of Xcq the ra dial distri- 
bution starts to resemble that fitted by llsraell (l2000t) . The off- 
sets of course will change if the normalization at /? ga i=8.5 kpc is 
adjusted . Note that there is a difference in the method u sed by 
IWilsonl (fl99l and lisrael (I2000T) to derive X C0 - IWilsonl dl995l) 
observed 12 CO in several giant molecular clouds within several 
nearby galaxies and compared their virial masses to the mass 
derived fro m thei r 12 CO luminosities using a Galactic conver- 
sion factor, llsraell (|2000) estimated the H2 mass of a sample of 
nearby galaxies by modelling their FIR continuum emission to 
derive their masses and related them to their 12 CO luminosities 
to derive Xco- The FIR continuum emission should trace the to- 
tal H2 content (including the CO-dark H2 component) while the 
method based on the virial mass derived from 12 CO mapping 
accounts for only the H2 mass traced by 12 CO 

The increase of Xqo with Galactocentric distance for the case 
when only CO and I3 CO are used to derive A^Fb.) is a result of 
the abundance gradient we applied in the calculation of N(H2). 
If we set the [CO]/[H2] abundance to be constant when con- 
verting JV(CO) to N(U 2 ) for Mask C02, and from A v to N{H 2 ) 
for MaskCOl, we obtain a nearly constant radial distribution 
of Xco- When the CO-dark H2 is included, however, the slope 
(seen in Figure l20li results from both the abundance and the ther- 
mal pressure gradients. If we also set the [C + ]/[H2] abundance to 
be constant we still obtain a distribution that increases with R ga \. 
Thus, the difference in the slopes in the distribution of .Xco is a 
result of the use of the thermal pressure gradient in the calcula- 
tion of the CO-dark H2 column density distribution. 

5.5. The surface density distribution and mass of the atomic 
and molecular gas in the Galactic plane 

In the top panel of Figure [21] we show the vertical surface 
density distribution of the total hydrogen gas in the Galactic 
plane. We also show the relative contributions from atomic 
and molecular gas. We calculated the surface densities assum- 
ing a Gaussian vertical distribution of the gas. For the molecu- 
lar gas we assume a constant FWHM disk thickness of 130pc 
(Section [3T2l i. In the case of the atomic gas, the disk thickness is 
observed to increase with Galactocent ric distance, and we use d 
the exponential relation presented by iKalberla & Kerpl (|2009), 
Zfwhm = 0.3e (s e"i _8 - 5 )/ 9 - 8 kpc. In all cases, the surface density 
includes a correction to account for the contribution from He. 
In the lower panel of Figure [21] we show the surface density 
distribution of the different contributing ISM components of the 
atomic and molecular gas studied here, the cold and warm neu- 
tral medium, the CO-dark H2 gas, and the H2 gas traced by 
I2 CO and 13 CO The surface density distribution of atomic and 
m olecular gas derived in our analysis is similar to that derived 
bv lScoville & Sanders! (1 19871) . but with the addition of the CO- 
dark H2 component, the distribution of molecular gas extends 
over a larger range of Galactocentric distances compared with 
that traced by 12 CO and 13 CO only. In Tabled we list the mass 
of each ISM component studied here. The CO-dark Kb gas ac- 
counts for about 30% of the total molecular gas mass. The WNM 
dominates the mass of atomic gas in the Galactic plane, account- 
ing for 70% of the total atomic mass. 



6. The FUV radiation field distribution in the Milky 
Way 

An important parameter governing the thermal balance and 
chemistry in interstellar clouds is the strength of the FUV field, 
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Fig. 21. {top panel) Radial distribution of the surface density of 
hydrogen in the Galactic plane. We also show the contribution 
of the atomic and molecular hydrogen gas. {lower panel) Radial 
distribution of the different ISM components, including the cold 
and warm neutral medium (atomic gas), the CO-dark H2 compo- 
nent, and the H2 gas that is traced by 12 CO and 13 CO. All surface 
densities have been corrected to account for the contribution of 
He. 



Xa (in units of the iDrainal 19781 field). The gas heating in cold, 
diffuse regions is dominated by the photoelectric effect resulting 
from the absorption of FUV photons by dust grains and PAHs. 
The FUV field also plays a critical role in the formation and de- 
struction of CO, determining the column density at which the 
C + /C°/CO transition takes place. It is therefore of interest to de- 
termine the typical strength of the FUV field to which the inter- 
stellar gas is exposed and how this FUV field varies from place 
to place in the plane of the Milky Way. 

PDR model calculations suggest that the [Cn]/ 12 CO inte- 
grated intensity ratio is prim arily a function of xo and the total H 
("hi + 2«h 2 ) volume density dWolfire et al.lll989tlKaufman et alJ 
1999). The reason for this dependence is that the 12 CO emission 
becomes optically thick quickly after a modest fraction of the 
gas-phase carbon is converted to CO. The [C n]/ 12 CO ratio is 
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Fig. 22. Comparison between the [Cn] and 12 CO emissivities 
corresponding to different Galactocentric rings, for gas associ- 
ated with Mask 3 (red circles) and Mask4 (black boxes). The 
straight lines represent fits to the data. Typical uncertainties are 
0.02 and 0.1 Kkms 1 kpc 1 , for [Cn] and 12 CO, respectively. 



thus determined by the column density of C + and the tempera- 
ture of the line-emitting region. Both quantities depend, in turn, 
on the strength of FUV field and H volume density. The solution 
for a given value of the [C n]/ I2 CO ratio is not unique, however, 
and several pairs of xo and «h can reproduce the observations. 
Therefore, we can only use the [C n]/ 12 CO ratio to determine a 
range of plausible values of the FUV radiation field and H vol- 
ume density that characterize the line-emitting region. 

In a preliminary study of [Cn], 12 CO, and 1 3 CO ve locity 
components on 16 GOTC+ LOSs, iPineda et all (1201 Obi) used 
PDR model calculations to constrain the strength of the FUV 
field to be mostly xo =1-10, over a wide range of H2 volume 
densities, thus suggesting that most of the observed [C 11] emis- 
sion associated with 12 CO and 13 CO is exposed to weak FUV 
radiation fields. A similar result, was obtained fro m PDR mod- 
eling o f the [Cn] emission observed by COBE bv lCubicket al.l 
( 2008), suggesting that most of the [C n] in our Galaxy originates 
in a clumpy medium exposed to an FUV field of xo — 60. 

With the complete GOTC+ survey, we can study the 
[C n]/ I2 CO ratio over the entire Galactic plane and determine 
how the strength of the FUV field is distributed in the Galaxy. 
We consider the azimuthally averaged emissivity of [C 11] cal- 
culated from spaxels where both 12 CO and 13 CO emission is 
detected (Mask 4), and where I2 CO is detected but 13 CO is not 
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Fig. 23. Radial distribution of the [Cn] to 12 CO emissivity ra- 
tio for Mask 3 (red circles) and Mask 4 (black boxes). We also 
show the radial distribution of the [C n]/ 12 CO ratio resulting 
from the combined emission from Mask 3 and 4 (blue triangles). 
The horizontal lines correspond to the average value for each 
Mask region between 3-10 kpc. The horizontal lines are labeled 
with a range of pairs of the FUV field and H2 volume densities 
(Xo, "h,X that repr oduce their cor r espon ding [C n]/ 12 CO ratio, 
as predicted by the lKaufman et alj d!999h PDR model. To facil- 
itate the comparison with the PDR model, the ratios were calcu- 
lated using [C11] and 12 CO emissivities expressed in units of erg 
s -1 cirT 2 sr~' kpc -1 . In these units, typical uncertainties in the ra- 
tios are 180, 33, and 35 for Mask 3, Mask 4, and the combined 
emission from both Mask regions, respectively. 



in volume densities because we do not expect a large difference 
in the spatial extent of such regions. We thus show the results of 
the PDR model for ranges in n^ that, based on the analysis in 
Section l5.2.2l might be appropriate. We considered the range of 
H 2 volume density 10 2 - 10 3 cirT 3 for Mask 3, and 10 3 - 10 4 crrT 3 
for Mask 4. For the combined emission from Mask 3 and Mask 4, 
we considered a wider range of possible H2 volume densities 
10 2 — 10 4 cm -3 . Note tha t the PDR model calculations presented 
bv lKaufman et alj i 19991) are in terms of the total H volume den- 
sity, including atomic and molecular hydrogen. For the compari- 
son between the observations and model predictions we assumed 
that the contribution of atomic hydrogen to the total H volume 
density is negligible and show the results of the PDR model in 
terms of the H2 volume density. 

We find that over the selected range of H2 volume densities, 
the [Cn]/ 12 CO ratio in both Mask 3 and Mask 4 suggest values 
of^o that are in the range 1-30, while we do not find any system- 
atic variation in xo with Galactocentric distance between 3 kpc 
and 10 kpc. The same range of^o is suggested by the [Cn]/ 12 CO 
ratio resulting from the combined emission from Mask 3 and 4, 
over a wider range of H2 volume densities. Our results therefore 
suggest that on average the [C 11] emission associated with 12 CO 
and/or I3 CO observed in the Galactic plane by the GOT C+ sur- 
vey emerges from gas exposed to modest FUV radiation fields, 
away from massive star-forming regions. 



7. Summary 

In this paper we have presented results of our study of the dis- 
tribution and properties of the different ISM components in the 
plane of the Milky Way using spatially and velocity resolved 
observations of the [C n] 158;um line. Our results can be summa- 
rized as follows: 



1. 



(Mask 3; See Figure[9]l. As discussed in Section[4] gas associated 
with Mask 3 might represent low-column density gas surround- 
ing the high-column density regions associated with Mask 4. In 
Figure [22] we show a comparison between the [Cn] and 12 CO 2. 
intensities for different Galactocentric rings, for gas associated 
with Mask 3 and 4. In both mask regions we see a good corre- 
lation between these two emissivities. A straight line fit results 3. 
in 7( 12 CO) = 7.8/([Cn]), for Mask 3, and 7( 12 CO) = 3.3/([Cn]) 
for Mask 4. 

In Figure [23] we show the radial distribution of the [C n] 
to 12 CO emissivity ratio for Masks 3 and 4. We also include 4. 
the radial distribution of the [C n]/ 12 CO ratio resulting from the 
combined emission from Masks 3 and 4. The ratio is computed 
considering emissivities in units of erg s -1 cirT 2 sr~' kpc" 1 . 
We restrict the range in Galactocentric distance to be between 
3 and 10 kpc because, for these regions, both [Cn] and 12 CO 
have significant emission in this range. As expected, given the 5. 
linear correlation between [C 11] and 12 CO, the ratio shows lit- 
tle variation between 3 and 10 kpc. We also include horizontal 
lines representing the average [C n]/ 12 CO ratio for each mask 
region. These horizontal lines are labeled with a range of pairs 
of the FUV field and H2 volume densities (xo, nn 2 ), that can 
reproduce the corresponding values of the [C n"| / 12 CO ratio, as 
predic ted by the PDR model calculations from iKaufman et alj 
(Il999l) . As mentioned before, Mask 3 might represent regions 6. 
with lower column densities compared to those in Mask 4. This 
column density difference might also translate into a difference 



We presented the first longitude-velocity maps of the [C 11] 
158/im line in the plane of the Milky Way. We found that 
the [C 11] emission is mostly associated with the spiral arms, 
tracing the envelopes of evolved clouds as well as clouds 
that are in the transition between the atomic and molecular 
phases. 

We found that most of the [C n] emission in the Galactic 
plane emerges from Galactocentric distances between 4 kpc 
and 11 kpc. 

We estimated that most of the observed [Cn] emission 
emerges from dense photon dominated regions (-47%), with 
smaller contributions from CO-dark H2 gas (-28%), cold 
atomic gas (-21%), and ionized gas (-4%). 
We used the [C 11] emission to separate the WNM and CNM 
contributions in the Galactic plane. We find that most of the 
atomic gas inside the Solar radius is in the form of cold neu- 
tral medium, while the warm neutral medium gas dominates 
the outer Galaxy. The average fraction of the total atomic gas 
that is CNM is -43%. 

We found that the warm and diffuse CO-dark H2 is dis- 
tributed over a larger range of Galactocentric distances (4- 
1 1 kpc) than the cold and dense H2 gas traced by 12 CO and 
13 CO (4-8 kpc). The fraction of the total H 2 gas that is CO- 
dark H2 increases with Galactocentric distance, ranging from 
-20% at 4 kpc up to -80% at 10 kpc. We estimate that the 
CO-dark H2 gas component accounts for -30% of the total 
molecular gas mass in the Milky Way. 
We estimated the radial distribution of the CO-to-H2 con- 
version factor. When N(H2) is derived from 12 CO and 
13 CO only, the CO-to-H2 conversion factor increases with 
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Galactocentric distance, influenced by the Galactic metallic - 
ity gradient. When the contribution of CO-dark H2 gas is in- 
cluded in NQ&i), the CO-to-H2 conversion factor increases 
with Galactocentric distance with an steeper slope, which is 
influenced by both metallicity and thermal pressure gradi- 
ents. 
7. We estimated the average strength of the FUV field mod- 
elling the [C n]/ 12 CO ratio. We find that most of the observed 
[C n] emission emerging from dense PDRs is associated with 
modest FUV fields in the xo - 1 - 30 range. 

Our results confirm that [C 11] is an important tool to char- 
acterize the different components of the ISM. Large-scale map- 
ping of the Milky Way and nearby galaxies with orbital and sub- 
orbital observatories are essential to further characterize the life- 
cycle of the interstellar medium, and will help us to build impor- 
tant templates for the interpretation of observations of [C 11] at 
high-redshift. 
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Appendix A: Azimuthally averaged line emissivities 

In the following we describe the procedure we used to calculate 
the radial distribution of the azimuthally averaged integrated in- 
tensity per unit length of a spectral line. Assuming a flat rotation 
curve and circular motions, the distance to the Galactic center R 
for a given data point with Galactic longitude I, latitude b, and 
local standard of rest (LSR) velocity Vlsr, is given by 



R = R a \l + 



Vlsr 



V Q sin(Z) cos(b) 



where R Q is the distance from the sun to the Galactic center and 
V Q is the orbital velocity of the sun with respect to the Galactic 
center. We adopt the IAU recommended values of R G - 8.5 kpc 
and V = 220 km s~'. We can invert Equation (1A.U to relate 
the observed LSR velocity of the emission to the Galactocentric 
distance through 



Vlsr = ^0 sin(0 cos 



«(*-.). 



(A.2) 



We divided the Galaxy in a set of concentric rings, with a 
central radial distance R and a radial extent AR - R out - R m . For 
a given ring, we used Equation ( !A.2t to calculate the correspond- 
ing interval limits in Vlsr- We then calculated the integrated 
intensity /(/, R) within this LSR velocity range. The integrated 
intensity depends on the path length crossed by a given line of 
sight, L(l, R out , R m ), and therefore we divided by this path length 
to calculate the integrated intensity per unit length (emissivity). 
The units of the emissivity are K km s~' kpc -1 . The value of L 
depends on the location of the Galaxy we are sampling. This de- 
pendency is illustrated in Figure lATI where we show a schematic 
of the different cases we considered for the calculation of L. 

In Case la we considered rings with outer radii inside the 
solar circle (i? ou t ^ Re) an d inner radius R[ n larger than the ra- 
dius of the tangent R tdn - \R Q sin(/)|. In this case a given ring is 
crossed twice by a given LOS. The path length for one of these 
crossings is given by, 



L = VCRout) 2 - ORtan) 2 - VORin) 2 - Wan) 2 . 



(A.3) 



In Case lb we considered rings with Ri n < R ma < R out . In 
this case a ring is crossed only once and the path length is given 

by, 



L = 2 VORout) 2 - Wan) 2 . 



(A.4) 



In Case 2 we considered LOSs that cross rings in the outer 
Galaxy with R\ n > R Q . These rings are crossed once and the path 
length is given by Equation (IA.31 >. 

Rings in Case 3a have R ln < R G < R oul and are observed in 
LOSs in the range 270°> / > 90°. In this case the path length is 
given by the distance to the point with R = R oa , 



L=R e COS(Z) + V(#out) 2 - Wan) 2 . 



(A.5) 



Finally for Case 3b, we considered rings with R m < R Q < 
Rout observed in LOSs in the range 270°< / < 90°. In this case 
we cross a ring twice. For the near side the path length is given 
by 



L=R Q cos(l)+ VCRin) 2 



Wan) 2 , 



(A.6) 



and for the far side is given by Equation (1A.3) . 

The final step is to average the emissivities over all samples 
of a given ring. This sum is done over all rings sampled by the 
GOT C+ survey N s , 



k(R) 



N, 



(A.7) 



We estimated the errors on the emissivity by propagating the 
errors of the individual determinations of the integrated intensity 
per unit length, AI/L, using, 



(A.l) AIq(R) 



^(AUIU) 2 



N, 



(A.8) 
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Case 3a 



Case 3b 




Case la 



Case 2 



Case lb 

Fig. A.l. Schematic of the Milky Way divided by a set of rings centered at the Galactic center. The two thick arrows originating 
from the Sun (shown as a circle) correspond to two lines-of-sight, one looking toward the inner Galaxy and the other toward the 
outer Galaxy. The dashed lines represent the distance to the Galactic center {R Q ) and the distance to the tangent point (R tdn ) for the 
line of sight that goes towards the inner Galaxy. The thin lines point to the different path lengths crossed by the lines-of-sight that 
correspond to the different cases considered in the text. 



In the top panel of Figure lAT2l we show the resulting ra- 
dial distribution of the [C n] emissivity for different values of 
AR - 0.5 kpc, 1 kpc, and 2kpc. We can see that for larger val- 
ues of AR the peak shifts to around 5 kpc, which is the midpoint 
where most of the [C n] emission observed in the GOT C+ sur- 
vey arises in the Galaxy. We adopt AR = 0.5 kpc to preserve 
the details of the [C n] emissivity distribution. The uncertainty 
in the [C n] emissivities for this selection of AR is typically 0.02 
K kms 1 kpc 1 . The bottom panel shows the number of sam- 
ples (N s ) for a given ring. As we can see, rings with radii below 
2 kpc are undersampled in our survey. Given that such rings are 
sampled by lines-of-sights with |/| < 3°, and that these LOSs 
are also known to have non-circular motions, we exclude these 
LOSs in our analysis. This exclusion has a minimal effect on 
emissivity distribution beyond 2 kpc. 



in Section gj] for the GOTC+ lines-of-sight with b 
and£ = ±1.0°. 



±0.5° 



Appendix B: Position-velocity maps 

In Figures IB. II to IB .41 we present the position-velocity maps 
of the [Cn], CO, and Hi emission for the GOTC+ lines-of- 
sight with b - ±0.5° and b - ±1.0°. Figure |B"31 shows the 
position-velocity maps for the different cloud types discussed 
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Fig. B.l. Position-velocity maps of the Milky Way in [C n] observed by GOT C+ for b - +0.5° 



:: 



100 

90 

80 , 

70 

60 

50 

40 

30 

20 

10 



co 

CD 
Q. 

E 




CO 



CD 



100 






2.5 



2 1 

CO 

CD 

a 
E 




1.5 



E 

CO 

CO 



100 







I I , I 

1 ■ 

i b = +0.5 deg 




- 




, Scutum-Crux — 







—»!-.'C~ 


, Sagitarius-Carina 

j :-'\t~\ Norma-Cygnus 






<:;::;:;;; 


:::::::„„ ■ ' 'l! J i Perseus 

'"a::;;. | . i 1 / 1 . 








: ' : ''"""•■■ ■".:,;,•• •'•'"'.ylll 1 JvfLH'* "■ • ' .'•*'' i 




- 


'••** 


x --;. : ; -.v;;^^%5'r':: ' 

; i iCi i ,| I li I.Uk^%^ 1 ' "'*^Z, ::::::::::: " -z, ::;:::;:;;;;; " 

. "-j-Li [--" : y|i ( .-- :.y . j -*.VC 


" 


- 




■ i !? : /)•' ii *•••>.. 


" 


- 


CO 


1 
1 1 1 


- 



0.5 



20 



18 



16 



14 £ 

CO 

12 CD 
Q. 

10I 

8 | 



6 ™ 



100 



24 



Pineda, Langer, Velusamy, and Goldsmith: A Herschel [C n] Galactic Plane Survey 



200 

150 

100 

50 



a o 

o 
cr -50 



-100 
-150 h 
-200 



HI 



-100 
200 

150 h 
100 
50 



I 
o 

CD 

or -50 

_l 

-100 
-150 
-200 

200 

150 
100 



-100 



50 - 



■B 

p 

CD 

rr -50 

03 



-100 



-150 - 



-200 



HU1 J ! : 



■ ■ : 



^S, 



I I 1 
.I 



'ill! ill? ! 



b = -0.5 deg 

Scutum-Crux 

■Sagitarius-Carina 

Norma-Cygnus 

Perseus 



■I 



r ■" 



4+ 

'ilMi 



„iM>t 



■.:..y« 

II 






.;,< 



jf l ' : 

Mil 



■I 



ffijhjrm 

"ffil'fT't'i ■t'k.iJ '■■■ 

••••fr^-j... .]..: ! i ! I r 



III! 



! 



-50 

Galactic Longitude [degrees] 



50 



-50 

Galactic Longitude [degrees] 



50 



b = -0.5 deg 

Scutum-Crux 

Sagitarius-Carina 

Norma-Cygnus 

Perseus 



- + .j. 



.iuM-I-'-'-Ij.,. , 

..;.:,■•-:':■;.;■■■... /.< y ■ ■. ... •• ■ i 

-. K«fel 

h — ■'v..-- I / ,f ••>-: 



CO 



-100 



-50 







50 



Galactic Longitude [degrees] 
Fig. B.2. Position-velocity maps of the Milky Way in [C n] observed by GOT C+ for b - -0.5° 
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Fig. B.3. Position-velocity maps of the Milky Way in [C n] observed by GOT C+ for b - + 1 .0° 
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Fig. B.4. Position-velocity maps of the Milky Way in [C n] observed by GOT C+ for b - - 1 .0° 
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Fig. B.5. Position velocity maps of the different cloud types defined in Section I4TI for b - +0.5° and +1.0°. The Mask {grey) 
represents velocity components with only H i detected, Mask 1 {white) are components with only H i and CO detected, Mask 2 
{blue) components with only H i and [C n], Mask 3 {black) components with H i, [C ii], and CO, and Mask 4 {Red) components with 
Hi, [Cn], 12 CO, and 13 CO. 
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